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Abstract 

A large number of treatments of the meson spectrum have been tried that consider mesons as 
quark - anti quark bound states. Recently, we used relativistic quantum "constraint" mechanics to 
introduce a fully covariant treatment defined by two coupled Dirac equations. For field-theoretic 
interactions, this procedure functions as a "quantum mechanical transform of Bethe-Salpeter equa- 
tion" . Here, we test its spectral fits against those provided by an assortment of models: Wisconsin 
model, Iowa State model, Brayshaw model, and the popular semi-relativistic treatment of God- 
frey and Isgur. We find that the fit provided by the two-body Dirac model for the entire meson 
spectrum competes with the best fits to partial spectra provided by the others and does so with 
the smallest number of interaction functions without additional cutoff parameters necessary to 
make other approaches numerically tractable. We discuss the distinguishing features of our model 
that may account for the relative overall success of its fits. Note especially that in our approach 
for QCD, the resulting pion mass and associated Goldstone behavior depend sensitively on the 
preservation of relativistic couplings that are crucial for its success when solved nonperturbatively 
for the analogous two-body bound-states of QED. 
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I. INTRODUCTION 



Over 50 years after the discovery of the first meson and over 25 years after the identifica- 
tion of its underlying quark degrees of freedom, the Strong-Interaction Bound-state problem 
remains unsolved. Perhaps eventually the full spectrum of mesonic and baryonic states will 
be calculated directly from Quantum Chromo dynamics via lattice gauge theory. This would 
require use of techniques that were unknown to the founding fathers of QED. For the present 
though, researchers have had to content themselves with attempts to extend bits and pieces 
of traditional QED bound-state treatments into the realm of QCD. Unfortunately, for those 
bound systems whose constituent kinetic or potential energies are comparable to constituent 
rest masses, nonrelativistic techniques are inadequate from the start. 

In the QED bound-state problem, weakness of the coupling permitted calculation through 
perturbation about the nonrelativistic quantum mechanics of the Schrodinger equation. Us- 
ing the equation adopted by Breit [1]- [3] (eventually justified by the Bethe-Salpeter equation 
[4]), one was faced with the fact that a nonperturbative numerical treatment of the Breit 
equation could not yield spectral results that agree to an appropriate order with a perturba- 
tive treatment of the semirelativistic form of that equation [4] -[9]. This form of the equation 
contained such terms as contact terms bred by the vector Darwin interaction that could be 
treated only perturbatively, spoiling the interpretation of the Breit equation as a bona-fide 
wave equation. Forays into the full relativistic structure defined by the Bethe-Salpctcr equa- 
tion turned up fundamental problems which fortunately could be sidestepped for QED due 
to the smallness of a. 

In the absence of definitive guidance from QED, in recent years researchers in QCD have 
felt free to jump off from any point that had proven historically useful in QED. Some have 
chosen to approach the spectrum using time-honored forms from the "relativistic correction 
structure" of atomic physics. Others have employed truncations of field-theoretic bound- 
state equations in hopes that the truncations do no violence to the dynamical structures 
or their relativistic transformation properties. A third set have broken away from QED 
by choosing to guess at "relativistic wave equations" as though such equations have no 
connection to field theory. 

Is there another way to attack this problem? Imagine that we could replace the 
Schrodinger equation by a many-body relativistic Schrodinger equation or improved Breit 
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equation that could be solved numerically. One would have to establish its validity by con- 
necting it to quantum field theory, and its utility by solving it for QCD. Of course such an 
approach would apply equally as well to QED and so would have to recapitulate the known 
results of QED. (These results might reemerge in unfamiliar forms since not originating in 
the usual expansion about the nonrclativistic limit.) 

Now, for the two-body bound-state problem, there is such an equation or rather a system 
of two coupled Dirac equations - for an interacting pair of relativistic spin one-half con- 
stituents. It turns out that for the two-body case, use of Dirac's constrained Hamiltonian 
mechanics [10]- [15] in a form appropriate for two spinning particles [18], [19] (pseudo-classical 
mechanics using Grassmann degrees of freedom [20], [21]) leads to a consistent relativistic 
quantum description. In the two-body case, one may explicitly construct the covariant 
Center of Momentum rest frame of the interacting system. In fact, the relativistic two- 
body problem may be written as an effective relativistic one-body problem [12], [22], [23]. 
The proper formulation of this relativistic scheme requires the successful treatment of the 
quantum ghost states (due to the presence of the "relative time") that first appeared in 
Nakanishi's work on the Bethe-Salpeter equation[24]. 

It might seem that although fully covariant and quantum mechanically legitimate, such 
an approach would merely give a sophisticated method for guessing relativistic wave equa- 
tions for systems of bound quarks. However, this method assumes its full power when 
combined with the field-theoretic machinery of the Bethe-Salpeter equation. When used 
with the kernel of the Bethe-Salpeter equation for QED, our approach combines weak- 
potential agreement with QED [25] with the nonperturbative structure of the field-theoretic 
eikonal approximation[26],[27]. The extra structure is automatically inherited from rela- 
tivistic classical[28] , [31] and quantum mechanics [2 7]. In QED our approach amounts to a 
"quantum- mechanical transform" [32], [33] of the Bethe Salpeter equation provided by two 
coupled Dirac equations whose fully covariant interactions are determined by QED in the 
Feynman Gauge[34],[25]. These "Two-Body Dirac Equations" are legitimate quantum wave 
equations that can be solved directly [35], [25] (without perturbation theory) whose numeri- 
cal or analytic solutions automatically agree with results generated by ordinary perturbative 
treatment. (In our opinion the importance of this agreement cannot be overemphasized. A 
common fault of most of the models we discuss in this paper is that they lack such agree- 
ment. But, if a numerical approach to a two body bound state formalism when specialized 
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to QED cannot reproduce the results given by its own perturbative treatment, how can 
one be certain that its apphcation to highly relativistic QCD bound states will not include 
spurious short range contributions.) 

Of course there is a fly in the ointment - but one to be expected on fundamental grounds. 
It turns out that the only separable interacting system as yet explicitly constructed in a 
canonical relativistic mechanics is the two body system. In practical terms, this means that 
we must confine the present treatment to the meson spectrum. So far, even the relativistic 
treatment of the three-body problem of QED in the constraint approach is unknown. No 
one has been able to produce three compatible separable Dirac equations which include not 
only mutual interactions but also necessary three body forces in closed form [13]. 

Although still considered unusual or unfamiliar by the bulk of bound-state researchers, the 
structures appearing in these equations may have been anticipated classically by J. L. Synge, 
the spin structures were introduced into QED (incorrectly) by Eddington and Gaunt [36], [37], 
and they have appeared in approximate forms appropriate for weak potentials in the works of 
Krapchev, Aneva, and Rizov[38] and of Pilkuhn[39] . Of greatest surprise but greatest value 
(to the authors), their perturbative weak-potential versions were uncovered in QED by J. 
Schwinger in his virial treatment of the positronium spectrum[40]. The associated relativistic 
mechanics transforms properly under spin-dependent generalizations of generators found by 
Pryce [42] Newton and Wigner[43]. The techniques for quantization go back to those found 
by Dirac[10], and apphed by Regge and Hanson to the relativistic top[44], by Nambu to the 
string, by Galvao and Teitelboim to the single spin one-half particle [21], and by Kalb and 
Van Alstine [ll]and by Todorov[12] to the pair of spinless particles. Their progenitors can 
be found in the bilocal field theories of Yukawa, Markov, Feynman and Gell-Mann as well as 
the myriad treatments of the relativistic oscillator beginning with the work of Schrodinger. 

In this paper, we will compare our latest results for the meson spectrum provided by 
Two-Body Dirac Equations with the corresponding results given by a representative sample 
of alternative methods. The present paper is not a detailed account of this method (already 
presented elsewhere - see Refs.[25][45]and references contained therein). Neither is it an 
attempt to conduct an even-handed or thorough review. Rather, its purpose is to show how 
such an organized scheme fares in the real world of calculation of a relativistic spectrum and 
to contrast its results with those produced by an array of approaches, each chosen on account 
of popularity or structural resemblances or differences with our approach. In this paper we 
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consider only approaches like ours that do not restrict themselves to the heavy mesons but 
attempt fits to the entire spectrum thus obtaining a more demanding comparison. (We do 
not consider here the myriad of partial spectral results for either the light or heavy mesons 
appearing in the recent literature). Where possible, we shall show how certain distinguishing 
features of the various approaches are responsible for success or failure of the resulting fits 
to the meson spectrum. Whether our equations ultimately prove correct or not, they have 
the virtue that they are explicitly numerically solvable without additional revisions, cutoffs, 
etc. unlike certain other approaches whose spectral consequences depend on ad hoc revision 
necessary for numerical solution. 

All of the treatments wc examine attempt to describe the interactions of QCD through 
the inclusion of spin-dependent interactions that in part first appeared as small corrections 
in atomic physics. All include relativistic kinematics for the constituents. One contributor 
to the use of such techniques [46] has even asserted that all of the alternative approaches that 
include relativistic kinematics are actually equivalent to the nonrelativistic quark model, so 
that the detailed relativistic structure of the interaction makes no difference to the bound 
state spectrum. However, as we shall see in a fully relativistic description with no extraneous 
parameters, the detailed relativistic interaction structure in fact determines the success or 
failure of a calculation of the full meson spectrum from a single equation. 

The order of the paper is as follows: 

First, in Section II, we review enough of the structures of our Two-Body Dirac Equations 
and their origins in relativistic constraint dynamics to make clear the equations that we 
are solving and the relativistic significances of the potential structures appearing in them. 
(Those readers who are already familiar with constraint dynamics might wish to go directly 
to the QCD applications of section III.) In Section III, we detail how we incorporate 
the interactions of QCD into our equations by constructing the relativistic version of the 
Adler-Piran static quark potential[47] that we use when we apply our equations to meson 
spectroscopy. In Section IV, we examine the numerical spectral results that are generated 
by this application of the Two-Body Dirac Equations. 

The feature of our approach that most distinguishes it from other more traditional two- 
body formalisms is its use of two coupled constituent equations (instead of one) containing 
two-body minimal substitution forms and related structures that incorporate the minimal 
interaction form of the original one-body Dirac equation. In Section V we rewite this form of 
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the Two Body Dirac Equations first as an equivalent one that incorporates the interactions 
through the kernel structures that appear in most older approaches and second as an equiva- 
lent form closely related to the Breit equation. In this section we examine how the relativistic 
interaction structures of the constraint approach lead even for QED to classifications of in- 
teraction terms that differ from the designations used in some of the other approaches. In 
Section VI, we examine an attempt to use the Salpeter Equation to treat the meson spec- 
trum: the Wisconsin Model of Gara, Durand, Durand, and Nickisch[48]. Although these 
authors try to keep relativistic structures, they ultimately employ weak-potential approxima- 
tions and structures obtained from perturbative QED in a non-perturbative equation (with 
no check to see that the procedure even makes nonpcrturbativc sense in QED itself). In 
Section VII we examine the Iowa State Model of Sommerer, Spence and Vary [49] which uses 
a new quasipotential approach for which, in contrast to the Wisconsin model, the authors 
check that the equation makes nonperturbative sense in QED at least for the positronium 
ground state. In Section VIII, we examine the Breit Equation Model of Brayshaw [50], 
which illustrates the sort of successful fit that one can still obtain when one is allowed to 
introduce confining interactions (into the Breit equation) through terms whose relativistic 
transformation properties are ambiguous. In Section IX, we look at the most popular treat- 
ment - the Semirelativistic Model of Godfrey and Isgur [51]. This model includes a different 
smearing and momentum dependent factor for each part of the various spin-dependent inter- 
actions. Although each interaction is introduced for apparently justifiable physical reasons, 
this approach breaks up (or spoils) the full relativistic spin structure that is the two-body 
counterpart of that of the one-body Dirac equation with its automatic relations among the 
various interaction terms. We examine this model to see how well our fully covariant set of 
two-body Dirac equations, employing only three potential parameters used in two different 
invariant interaction functions, can do versus Godfrey's and Isgur's semirelativistic equation 
with relativistic kinematics and pieces of relativistic dynamical corrections (introduced in 
a patchwork manner with ten potential parameters used in six different interaction func- 
tions), when required to fit the whole meson spectrum (including the light-quark mesons). 
Finally, in Section X, we conclude the paper by reviewing some of the features of the con- 
straint approach that played important roles in the relative success of its fit to the meson 
spectrum. We then use apparent successes of recent fits produced by the ordinary nonrela- 
tivistic quark model to point out dangers inherent in judging rival formalisms on the basis of 
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fits to portions of the spectrum. At the end of the paper, we supply sets of tables for spectral 
comparisons and appendices detailing the radial form of our Two-Body Dirac equations that 
we use for our spectral calculations, and the numerical procedure that we use to construct 
meson wave functions. We also include a table summarizing the important features of the 
various methods that we compare in this paper. 

II. THE TWO-BODY DIRAC EQUATIONS OF CONSTRAINT DYNAMICS 

In order to treat a single relativistic spin-one-half particle, Dirac originally constructed 
a quantum wave equation from a first-order wave operator that is the matrix square-root 
of the corresponding Klein-Gordon operator [52]. Our method extends his construction to 
the system of two interacting relativistic spin-one-half particles with quantum dynamics 
governed by a pair of compatible Dirac equations on a single 16-componcnt wave function. 
For an extensive review of this approach, see Refs.[23, 25, 45] and works cited therein. For 
the reader unfamiliar with this approach, we present a brief review. 

About 27 years ago, the relativistic constraint approach first successfully yielded a co- 
variant yet canonical formulation of the relativistic two-body problem for two interacting 
spinless classical particles. It accomplished this by covariantly controlling the troublesome 
relative time and relative energy, thereby reducing the number of degrees of freedom of the 
relativistic two-body problem to that of the corresponding nonrelativistic problem[ll]-[15]. 
In this method, the reduction takes place through the enforcement of a generalized mass shell 
constraint for each of the two interacting spinless particles: pf + mf + ^i 0. Mathematical 
consistency then requires that the two constraints be "compatible" in the sense that they 
be conserved by a covariant systcm-Hamiltonian. Upon quantization, the quantum version 
of this "compatibility condition" becomes the requirement that the quantum versions of the 
constraints (two separate Klein-Gordon equations on the same wave function for spinless 
particles) possess a commutator that vanishes when applied to the wave-function. In 1982, 
the authors of this paper used a supersymmetric classical formulation of the single-particle 
Dirac equation due to Galvao and Teitelboim to successfully extend this construction to 
the "pseudoclassical" mechanics of two spin-one-half particles [18, 21]. Upon quantization, 
this scheme produces a consistent relativistic quantum mechanics for a pair of interacting 
fermions governed by two coupled Dirac equations. 
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When specialized to the case of two relativistic spin-one-half particles interacting through 
four-vector and scalar potentials, the two compatible 16-component Dirac equations [23, 25, 
45] take the form 

Si^l^ = 751 (71 • (Pi - ^i) + ^1 + Si)jIj = (2.1a) 

S2IP = 752(72 ■ (P2 - A2) +1712 + S2)lp = 0, (2.1b) 

in terms of Si operators that in the free-particle limit become operator square roots of the 
Klein-Gordon operator. 

The relativistic four- vector potentials A'^ and scalar potentials Si are effective constituent 
potentials that in either limit — > 00 go over to the ordinary external vector and scalar 
potentials of the light-particle's one-body Dirac equation. Note that the four-vector interac- 
tions enter through "minimal substitutions" inherited (along with the accompanying gauge 
structure) from the corresponding classical field theory[22, 28, 31]. The covariant spin- 
dependent terms in the constituent vector and scalar potentials (see Eq.(2.10 and Eq.(2.11) 
below) are recoil terms whose forms are nonperturbative consequences of the compatibility 
condition 

[Si,S2]^ = 0. (2.2) 

This condition also requires that the potentials depend on the space-like interparticle sepa- 
ration only through the combination 

= + W'')(a;i-X2). (2.3) 

with no dependence on the relative time in the cm. frame. This separation variable is 
orthogonal to the total four-momentum 

PM^pM+pM. _p2^^2 ^2.4) 

P is the time-like unit vector 

= P7w. (2.5) 
The accompanying relative four-momentum canonically conjugate to x± is 

p'^ = {€2P2 — eiP2)/w;wh.eie €1 + €2 = w, €1 — €2 = {ml — m\)/w (2.6) 

in which w is the total cm. energy. The e^'s are the invariant cm. energies of each of the 
(interacting) particles [29]. 
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The wave operators in Eqs.(2.1a,2.1b) operate on a single 16-component spinor which we 
decompose as 

= (2.7) 

^3 

in which the are four-component spinors. 

Once we have ensured that the compatibihty condition is satisfied, Eqs.(2.1a,2.1b) provide 
a consistent quantum description of the relativistic two-body system incorporating several 
important properties [23, 25, 45] . They are manifestly covariant. They reduce to the 
ordinary one body Dirac equation in the limit in which either of the particles becomes 
infinitely heavy. They can be combined to give [25, 53] coupled second-order Schrodinger-like 
equations (Pauli-forms) for the sixteen component Dirac spinors. In the center of momentum 
(cm.) system, for the vector and scalar interactions of Eq.(2.10) and Eq.(2.11) below, these 
equations resemble ordinary Schrodinger equations with interactions that include central- 
potential, Darwin, spin-orbit, spin-spin, and tensor terms. These customary terms are 
accompanied by others that provide important additional couplings between the upper-upper 
{ipi) and lower- lower ('^4) four component spinor portions of the full sixteen component Dirac 
spinor. The interactions are completely local but depend explicitly on the total energy w in 
the cm. frame. In this paper we use a recently developed rearrangement of these equations 
[53] (similar to that first presented in [54]) that provides us with ones simpler to solve 
but physically equivalent The resulting local Schrodinger-like equation depending on the 
four- component spinor (p^ = ipi + ■?/'4 takes the general cm. form 

(_V2 + $(r, o-i, a„ w))(t>+ = h\w)4>+. (2.8) 

with no coupling to other four component spinors. The explicit version of the potential $ in 
Eq. (2.8) that results from the rearrangement has a structure that produces couplings between 
the spin components of 0+ that are no more complicated than those of its nonrelativistic 
counterpart - with the customary spin-spin, spin-orbit, non-central tensor or spin-orbit 
difference terms appearing. We have checked that both the simpler form Eq.(2.8) and 
the equivalent coupled forms give the same numerical spectral results when tested for QED 
bound states as in [25] and when tested for our new QCD spectral results appearing in this 
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paper. (This provides an important cross check on our numerical calculation of the meson 
spectra). Eq.(2.8) is accompanied by similar equations for 0_ = ■^i — ■04 and x± = '02 ± '03- 
Once Eq.(2.8) is solved, one can use Eq.(2.1a,2.1b) to determine 0_ and x±- Because of 
the decoupling it is not necessary to determine 0_ and x± to solve the eigenvalue equation 
(2.8). However, the detailed form of $ for 0+ results from their elimination through the 
Pauli reduction procedure. In these equations, the usual invariant 

b'^{w) = (w^ - 2w^{ml + ml) + {mj - m^)2)/4w^ (2.9) 

plays the role of energy eigenvalue. This invariant is the cm. value of the square of the 
relative momentum expressed as a function of the invariant total cm. energy w. 

Note that in the limit in which one of the particles becomes very heavy, this Schrodinger- 
like equation turns into the one obtained by eliminating the lower component of the ordinary 
one-body Dirac equation in terms of the other component. 

The vector potentials appearing in Eqs.(2.1a,2.1b) depend on three invariant functions 
El, E2, and G that define time-hke vector interactions (proportional to P) and space-hke 
vector interactions (orthogonal to P, with = d/dx^) [23, 25] 

C riF' i 

At = ((61 - El) - i-72 • (-^ + aG)72 ■P)P^+{1- G)p^ - -dG ■ 7272^ 

G r)E 7 
Ai^ = {{e, - E,) + z-71 • (-^ + aG)7i ■P)P'-{1- G)p^ + -dG ■ 717^, (2.10) 

while the scalar potentials Si depend on G and two additional invariant functions Mi and 
M2 

Si-Mi-mi--G^,— 

S2 = M2-m2+'-G^i-^. (2.11) 

Note that the terms in 2.10 and 2.11 which are explicitly spin-dependent through the gamma 
matrices are essential in order to satisfy the compatibility condition 2.2. Later on, when 
the equation is reduced to second-order "Pauli- form" , yet other spin dependences eventually 
arise from gamma matrix terms (that, when squared, lose their gamma matrix dependence). 
These are typical of what occurs in the reduction of the one-body Dirac equation to the 
"Pauli form" . The gamma matrices also give rise to spin independent terms in the Pauli- 
forms. These terms emerge in a manner similar to the above two sources of spin dependent 
terms in the Pauli-form of the equations. 
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In the case in which the space-hke and time-hke vectors are not independent but combine 
into electromagnetic-hke four-vectors, the constituent vector interactions appear in a more 
compact form 

= (e: - ^^^^^ + + (1 - G)p" - . T27? 

= (e, - + ^^j^"* - (1 - G),f + \dG ■ T.7r (2.12) 

In that case Ei, E2 and G are related to each other [22, 23] {dEi/ E2 = —d\ogG) and for our 
QCD appUcations (as well as for QED) are functions of only one invariant function A{r) in 
which r is the invariant 

r = J^^. (2.13) 



They take the forms 



El{A)^G\e,-Af, 

El{A) = G\e2 - A)\ (2.14) 



in which 



= 7 ^T^. (2.15) 

(1 - 2A/w) ^ ' 

In the forms of these equations used below, Todorov's collective energy variable 

= (w^ -ml- ml) /2w, (2. 16) 

will eventually appear. 

In general Mi and M2 are related to each other[18, 23] and for QCD apphcations are 
functions of two invariant functions A{r) and S{r) appearing in the forms: 

Ml{A, S)^m\^ G\2m^S + S^) 

M^{A,S) = ml + G\2m^S + S^), (2.17) 

in which 

m^ — mim2/w. (2-18) 

(In these equations, m^j and are the relativistic reduced mass and energy of the fictitious 
particle of relative motion introduced by Todorov [12, 26], which satisfy the effective one- 
body Einstein condition 

el-ml = b\w). (2.19) 
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In the limit in which one of the particles becomes infinitely heavy, mw and reduce to the 
mass and energy of the lighter particle.) The invariant function S{r) is primarily responsible 
for the constituent scalar potentials since Si — ii S{r) = , while A{r) contributes to 
the Si (if S{r) ^ 0) as well as to the vector potentials A^. Originally, we derived the 
general forms of Eqs. (2. 17,2. 14, 2. 15) for the scalar and vector potentials using classical 
field theoretic arguments [28, 31] (see also [18, 26]). Surprisingly, the resulting forms for 
the mass and energy potential functions Mj, G and Ei automatically embody collective 
minimal substitution rules for the spin-independent parts of the Schrodinger-like forms of 
the equations. Classically those forms turn out to be modifications of the Einstein condition 
for the free effective particle of relative motion 



For the vector interaction they automatically generate the replacement of by — ^ and 
for the scalar interaction the replacement of rriyj by + S. The part of Eq.(2.8) that 
results from the vector and scalar interactions then takes the form 



Now, we originally found these forms starting from relativistic classical field theory. The 
deceptively simple form of Eq.(2.21) in fact incorporates retarded and advanced effects 
through its dependnce on the cm. energy w. On the other hand, recently Jallouli and 
Sazdjian [27] obtained Eqs. (2. 14) and (2.17) in quantum field theory after performing a 
necessarily laborious Eikonal summation to all orders of ladder and cross ladder diagrams 
together with all constraint diagrams (Lippmann-Schwinger like iterations of the simple 
Born diagram) [55]. Thus, the structure first discovered simply in the correspondence limit 
has now been verified through direct but difficult derivation from perturbative quantum field 
theory. 

These equations contain an important hidden hyperbolic structure (which we could have 
used to introdiicc the interactions in the first place). To employ it wc introduce two inde- 
pendent invariant functions ^(xj,) and ^(xj,), in terms of which the invariant functions of 
Eqs.(2.10,2.11) take the forms: 



2,2 2 



(2.20) 



(p2 + 2m^S + S^ + 2e^A - A^)(p+ = h^. 



(2.21) 



Ml = mi cosh L + m2 sinh L 



M2 — 1712 cosh L +mi sinh L 



(2.22) 
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El — ei cosh 0—62 sinh Q 



E2 = 62 cosh g — 61 sinh Q 



(2.23) 



exp^. 



(2.24) 



In terms of Q and the constituent momenta pi and p2 , the individual four- vector potentials 
of Eq.(2.14) take the suggestive forms 

^1 = [1 - cosh(^)]pi + sinh(^)p2 - ^zidexpG ■ 72)72 



Eqs.(2.22), (2.23) and (2.25) together display a further consequence of the compatibility 
condition, a kind of relativistic version of Newton's third law in the sense that the two sets 
of constituent scalar and vector potentials are each given in terms of just one invariant 
function, S and A respectively. 

In terms of these functions the coupled two-body Dirac equations then take the form 




(2.25) 




in which 




depending on gamma matrices with block forms 
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Si = i^i(3au (2.34) 

As described in Appendix A, a procedure analogous to the Pauli reduction procedure of 
the one-body Dirac equation case yields 



r ^ 2E2Mi + E2Mi^ 2r E2M1 + E2Mi^ ' 

^2 2r 2^ ' 2^ 'SjMi + BjAf, 

^ ^ ^2 2r 2 ^ ' E2M1 + E2M1' 
+ 2 + (^1 • ra2 • p + ^2 • mi . p) + -{L - ^)^^^^^— ^ ]0+ (2.35) 

= 6'(^/;)0+ 

Eq.(2.35) is the coupled four-component Schrodinger-like form of our equations that we 
use for our quark model bound state calculations for the mesons in the present paper. 
It can be solved nonperturbatively not only for quark model calculations but also for QED 
calculations since in that case every term is quantum-mechanically well defined (less singular 
than -l/4r^). 

Prom this equation we obtain two coupled radial Schrodinger-like equations in the general 
case. But for j' = or spin singlet states these equations reduce to uncoupled equations. The 
extra component for the general case arises from orbital angular momentum mixing or spin 
mixing, the latter absent for equal mass states. The detailed radial forms of these equations 
are given in Appendix A. Por the case of QED { S = 0, A = —a/r), we have solved these 
coupled Schrodinger-like equations numerically obtaining results that are explicitly accurate 
through order cc^ (with errors on the order of q;^)[25]. We have even obtained analytic 
solutions to the full system of coupled 16 component Dirac equations in the important 
case of spin-singlet positronium [35]. Por both numerical and analytic solution, the results 
agree with those produced by perturbative treatment of these equations and with standard 
spectral results 
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III. MESON SPECTROSCOPY 



We use the constraint Eq.(2.35) to construct a relativistic naive quark model by choosing 
the two invariant functions Q and L or equivalently A and S to incorporate a version of 
the static quark potential originally obtained from QCD by Adler and Piran [47] through 
a nonlinear effective action model for heavy quark statics. They used the renormalization 
group approximation to obtain both total flux confinement and a linear static potential 
at large distances. Their model uses nonlinear electrostatics with displacment and electric 
fields related through a nonlinear constitutive equation with the effective dielectric constant 
given by a leading log log model which fixes all parameters in their model apart from a mass 
scale A. Their static potential contains an infinite additve constant which in turn results in 
the inclusion of an unknown constant Uq in the final form of their potential (hereafter called 
VAp{r)). We insert into Eq.(2.35) invariants A and S with forms determined so that the 
sum A + S appearing as the potential in the nonrelativistic limit of our equations becomes 
the Adler-Piran nonrelativistic QQ potential (which depends on two parameters A and Uq) 
plus the Coulomb interaction between the quark and antiquark. That is, 

VAp{r) + V^i = A(C/(Ar) +Uo) + ^^A+S, (3.1) 

r 

As determined by Adler and Piran, the short and long distance behaviors of U (Ar) generate 
known lattice and continuum results through the explicit appearance of an effective run- 
ning coupling constant in coordinate space. That is, the Adler-Piran potential incorporates 
asymptotic freedom through 

AC/(Ar « 1) ~ l/(rlnAr), (3.2) 

and linear confinement through 

AU{Ar >> 1) - AV. (3.3) 

The long distance ( = Ar > 2) behavior of the static potential VAp(r) is given explicitly by 

A{cix + C2 ln(a;) + -^ + — + 05) (3.4) 

in which x — Ar while the coefficients q are given by the Adler-Piran leading log-log model 
[47] . In addition to obtaining these analytic forms for short and long distances they converted 
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the numerically obtained values of the potential at intermediate distances to a compact 
analytic expression. The nonrelativistic analysis used by Adler and Piran, however, does 
not determine the relativistic transformation properties of the potential. How this potential 
is apportioned between vector and scalar is therefore somewhat, although not completely, 
arbitrary. In earlier work [34] we divided the potential in the following way among three 
relativistic invariants V(r), S, and ^.(In our former construction, the additional invariant V 
was responsible for a possible independent time-like vector interaction.) 



S = rj{A{cix + C2 ln(x) + ^ + cg + C/q), 

V = (1 - 7]){A{ciX + C2ln(a;) + ^ + Cg + t/o), 

A^Va-S-V, (3.5) 

in which r] — \. That is, we assumed that (with the exception of the Coulomb-like term 
{ci/x)) the long distance part was equally divided between scalar and a proposed time- 
like vector. In the present paper we drop the time-like vector for reasons detailed below 
and assume instead that the scalar interaction is solely responsible for the long distance 
confining terms (77 = 1). The attractive (04 = —0.58) QCD Coulomb- like portion (not to 
be confused with the electrostatic Koui) is assigned completely to the "electromagnetic-like" 
part A. That is, the constant portion of the running coupling constant corresponding to the 
exchange diagram is expected to be electromagnetic-like. 

Elsewhere we have treated another model explicitly containing these features: the 
Richardson potential. Its momentum space form 

V{q) - l/q^n n(l + qVA^) (3.6) 

interpolates in a simple way between asymptotic freedom V{q) ~ l/q^Zn(q^/A^) and linear 
confinement V^(q) ~ l/q'^- Even though the Richardson potential is not tied to any field 

theoretic base in the intermediate region (unlike the Adler-Piran potential) and docs not 
give as good fits to the data, it does provide a convenient form for displaying our points 
about the static quark potential. The Richardson radial form is 

V{r) = 87rAV/27 - 87r/(Ar)/(27r) (3.7) 
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For r — > 0, /(Ar) — > — l/ln(Ar), while for r — > oo, /(Ar) — > 1. Thus, in this model, if 
the confining part of the potential is a world scalar, then in the large r limit the remain- 
ing portion, regarded as an electromagnetic-like interaction corresponding to our invariant 
function A{r), would be an attractive 1/r potential with a coupling constant on the order 
of 1. This is in reasonable agreement with the Adler model which also has an attractive 
1/r part. Support for the assumption that the C4 term belongs only to A also arises from 
phenomenological considerations. We find that attempts to assign the C4 term to the scalar 
potential have a drastic effect on the spin-spin and spin-orbit splittings. In fact, using this 
term in S through Eqs.(2.17) generates spin-spin and spin-orbit splittings that are much too 
small. 

In our previous work, we divided the confining part equally between scalar and time-like 
vector so that the spin-orbit multiplets would not be inverted. This was done in order to 
obtain from our model the ao(980) meson which was then considered as the prime candidate 
for the relativistic counterpart of the ^Pq meson. However, recent analysis indicates that 
that meson may be instead a meson-meson or four quark bound state (see however, [59], 
which even interprets the ao(980) meson as part of a new scalar {^Sq) meson qq multiplet 
outside of the usual quark model) while a meson with mass of 1450 MeV may be the correct 
candidate for the quark model state [57]. Interpretation of this other state as the ^Pq meson 
would in fact require a partial inversion of the spin-orbit triplet (from what one would 
expect based on the positronium analog). This partial inversion is consistent with the ^Pq 
candidate for the us system also appearing in a position that partially inverts the spin-orbit 
splitting. Since the sole purpose of including V in our previous treatment was to prevent 
the inversion, we exclude it from our present treatment. In our older treatment [34], we 
neglected the tensor coupling, unequal mass spin-orbit difference couplings, and the u — d 
quark mass differences. In the present treatment, we treat the entire interaction present 
in our equations, thereby keeping each of these effects. In our former treatment we also 
performed a decoupling between the upper-upper and lower-lower components of the wave 
functions for spin-triplet states which turned out to be defective but which we subsequently 
corrected in our numerical test of our formalism for QED [25]. The corrected decoupling 
(appearing in Eq.(2.35)) is included in the new meson calculations appearing in this paper. 

In the present investigation, we compute the best fit meson spectrum for the following 
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apportionment of the Adler-Piran potential: 

A = exp(-/?Ar)[V^AP - ^)] + ^ + (3.8a) 

r r r 

S^Vap+^-^-A^{Vap- -)(1 - exp(-/?Ar)) (3.8b) 
r r 

In order to covariantly incorporate the Adler-Piran potential into our equations, we treat 
the short distance portion as purely electromagnetic-like (in the sense of the transformation 
properties of the potential). Through the additional parameter /3, the exponential factor 
gradually turns off the electromagnetic-like contribution to the potential at long distance 
except for the 1/r portion mentioned above, while the scalar portion gradually turns on, 
becoming fully responsible for the linear confining and subdominant terms at long distance. 
Altogether our two invariant potential functions depend on three parameters: A, C/q, and (5. 

When inserted into the constraint equations, S and A become relativistic invariant func- 
tions of the invariant separation r = \/x\^ . The covariant structures of the constraint 
formalism then embellish the central static potential with accompanying spin-dependent 
and recoil terms. 

In general applications of these two-body Dirac equations one must ensure that the 
values assumed by A and S always result in real interaction functions Ei,Mi, and G while 
preserving the correct heavy particle limits. In particular a large repulsive A will give an 
imaginary G while a large attractive S would lead in the limit when one particle becomes 
heavy to an incorrect form of the one-body Dirac equation ( for m2 oo the interaction 
mass potential function Mi \mi + S\ instead of mi + S). In the calculations contained 
in the present paper, the best fit parameters turn out to be such that A always remains 
attractive while S always remains replusive so we need not make any modifications. Such 
problems do arise in the nucleon-nucleon scattering problem. See [60] for a discussion of 
these problems and their resolution. 

IV. NUMERICAL SPECTRAL RESULTS 

A. Tabulation and Discussion of Computed Meson Spectra 

We now use our formahsm as embodied in Eqs. (2.35), and (3.8a,3.8b) to calculate the 
full meson spectrum including the light-quark mesons. (As a check on these calculations 
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we have also used the older forms derived in [25]). Note that the nonrelativistic quark 
model when used in conjunction with realistic QCD potentials such as Richardson's po- 
tential or the Adler-Piran potential fails for light mesons since the ordinary nonrelativistic 
Schrodinger equation's lack of relativistic kinematics leads to increasing meson masses as 
the quark masses drop below a certain point [64], thereby spoiling proper treatment of the 
pion, as well as other states. Here, we shall see how our relativistic equations remedy this 
situation. In addition to including the proper relativistic kinematics, our equations also 
contain energy dependence in the dynamical quasipotential. Mathematically, this feature 
turns our equations into wave equations that depend nonlinear ly on the eigenvalue. Their 
solution, which we have treated in detail elsewhere (see [34, 61]), requires an efficient iter- 
ation scheme to supplement our algorithm for the eigenvalue b'^{w) when our equations are 
written as coupled Schrodinger-like forms. 

We display our results in Table I at the end of the paper. In the first two columns of each of 
the tables we list quantum numbers and experimental rest mass values for 89 known mesons. 
We include all well known and plausible candidates hsted in the standard reference ([57]). We 
omit only those mesons with substantial flavor mixing. In the tables, the quantum numbers 
listed are those of the 0+ part of the sixteen-component wave function. To generate the fits, 
in addition to the the quark masses we employ the parameters A, C/q and We merely insert 
the static Adler-Piran potential into our relativistic wave equations just as we have inserted 
the Coulomb potential A — —a/r to obtain the results of QED[25, 35]. Note especially 
that we use a single ^{A, S) for all quark mass ratios - hence a single structure for all the 
QQ, Qq, and qq mesons in a single overall fit. In the third column in Table I we present the 
results for the model defined by Eqs.(3.8a,3.8b). The entire confining part of the potential 
in this model transforms as a world scalar. In our equations, this structure leads to linear 
confinement at long distances and quadratic confinement at extremely long distances (where 
the quadratic contribution outweighs the linear term 2m^S'). At distances at which 
exp(— /3Ar) << 1, the corresponding spin-orbit, Thomas, and Darwin terms are dominated 
by the scalar interaction, while at short distances (exp(— /3Ar) ~ 1) the electromagnetic-hke 
portion of the interaction gives the dominant contribution to the fine structure. Furthermore 
because the signs of each of the spin-orbit and Darwin terms in the Pauli-form of our Dirac 
equations are opposite for the scalar and vector interactions, the spin-orbit contributions of 
those parts of the interaction produce opposite effects with degrees of cancellation depending 
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on the size of the quarkonium atom. 

We obtain the meson masses given in column three as the result of a least squares fit 
using the known experimental errors, an assumed calculational error of 1.0 MeV, and an 
independent error conservatively taken to be 5% of the total width of the meson. We 
employ the calculational error not to represent the uncertainty of our algorithm but instead 
to prevent the mesons that are stable with respect to the strong interaction from being 
weighted too heavily. Our is per datum (89) minus parameters (8). 

The resulting best fit turns out to have quark masses — 4.877, — 1.507, = 
0.253, vriu — 0.0547, = .0580 GeV , along with potential parameters A = 0.216, AC/q — 
1.865 GeV and inverse distance parameter (5 = 1.936. This value of /? implies that (in the 
best fit) as the quark separation increases, our apportioned Adler-Piran potential switches 
from primarily vector to scalar at about 0.5 fermi. This shift is a relativistic effect since 
the effective nonrelativistic hmit of the potential {A + S) exhibits no such shift (i.e., by 
construction P drops out). 

In Table I, the numbers given in parentheses to the right of the experimental meson 
masses are experimental errors in MeV. The numbers given in parentheses to the right of 
the predicted meson masses are the contribution of that meson's calculation to the total 
of 101 . 

The 17 mesons that contain a b (or b) quark contribute a total of about 5.4 to the x^, at 
an average of about 0.3 each. This is the lowest contribution of those given by any family. 
Since the Adler-Piran potential was originally derived for static quarks, one should not be 
surprised to find that most of the best fit mesons are members of the least relativistic of the 
meson families. Note, however, that five of the best fit mesons of this type contain highly 
relativistic u and s quarks (for which our equation reduces essentially to the one-body Dirac 
equation for the light quark). 

The 24 mesons that contain a c (or c) quark contribute a total of about 50.7 to the at 
an average of about 2.2 each. This is the highest contribution of those given by any family. 
A significant part of this contribution is due to the ip meson mass being about 32 MeV off 
its experimental value. Another part of the contribution is due to fact that the mass of the 
high orbital excitation of the D* tensor meson is 80 MeV below its experimental value. In 
addition, the high orbital excitation of the D* is 60 MeV low. 

The 24 mesons that contain an s (or s) quark contribute a total of about 46.3 to the 
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at an average of about 1.3 each, less than that for the c-quark mesons. This is important 
because the s quarks are hghter than the c quarks. Part of the reason for this unexpected 
effect is that our fitting procedure accounts for the fact that our meson model ignores the 
level shifts (due to the instability of many of the mesons that contain an s quark) through 
the introduction of a theoretical error on the order of 5% of the width of the unstable mesons. 

The 36 mesons that contain a u (or u) quark contribute a total of about 54.6 to the at 
an average of about 1.5 each while the 16 mesons on our list that contain a d (or d) quark 
contribute a total of about 18.6 to the at an average of about 1.2 each. 

The worst fits produced by our model are those to the and the D* and D* high orbital 
excitations. Although two of these mesons contain the light u and d, quarks, in our fit the 
more relativistic bound states arc not in general fit less well. In fact, the tt, K, D and p 
mesons are fit better than these two excited D* and D* mesons. 

We see that over all, the two-body Dirac equations together with the relativistic version 
of the Adler-Piran potential account very well for the meson spectrum over the entire range 
of relativistic motions, using just the two parametric functions A and ^S". 

We now examine another important feature of our method: the goodness with which our 
equations account for spin-dependent effects (both fine- and hypcrfinc- splittings). Table 
I shows that the best fit versus experimental, ground state singlet-triplet splittings for the 
bu, bs, cc, cu, cd, cs, su, sd, ud systems are 48 vs 46, 59 vs 47, 151 vs 117, 134 vs 142, 
132 vs 142, 147 vs 144, 418 vs 398, 418 vs 398, and finally 648 vs 627 MeV. We obtain a 
uniformly good fit for all hyperfine ground state splittings except for the rjc — ip system. 
The problem with the fit for that system of mesons occurs because the D* and D* 
states arc significantly low while the ip is significantly high. Furthermore, the singlet and 
triplet P states are uniformly low. An attempt to lower the c quark mass by correcting the 
■0 mass while raising the charmonium and D*, D* P state masses would require raising the 
c quark mass. Reducing one discrepancy would worsen the other. Below, we will uncover 
what we believe is the primary cause for this discrepancy as we examine other aspects of 
the spectrum. 

For the spin-orbit splittings we obtain values for the R ratios {^P2 — ^ Pi)/(^Pi Pq)) 
of 0.71,0.67,0.42,-0.19,-0.58,-3.35 for the two hh triplets, and the cc, ss,us,ud spin triplets 
compared to the experimental ratios of 0.66,0.61,0.48,0.09,-0.97,-0.4. This fit ranges from 
very good in the case of the light T multiplet to miserably bad for the two hghtest multiplets. 
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Prom the experimental point of view some of the problem may be blamed on the uncertain 
status of the ^Pq light quark meson bound states and the spin-mixing in the case of the K* 
multiplet. Prom the theoretical point of view, the lack of any mechanism in our model to 
account for the effects of decay rates on level shifts undoubtedly has an effect. Another likely 
cause is that as we proceed from the heavy mesons to the light ones, the radial size of the 
meson grows so that the long distance interactions, in which the scalar interaction becomes 
dominant, play a more important role. The spin-orbit terms due to scalar interactions are 
opposite in sign to and tend (at long distance) to dominate the spin-orbit terms due to 
vector interactions. This results in partial to full multiplet inversions as we proceed from 
the ss to the ud mesons. This inversion mechanism may also be responsible for the problems 
of the two orbitally excited D* and D* mesons described above. It may be responsible as 
well for the problem of the singlet and triplet P states since the scalar interaction tends to 
offset the dominant shorter range vector interaction, at least slightly. 

We also examine the effect of the hyperfine structure of our equations on the splitting 
between the ^Pi and weighted sum [5{^P2) + 3(^-Pi) + l(^Po)]/9 of bound states. We obtain 
pairs of values equal to 3. 520,3.520;!. 408, 1.432; 1.392,1416 for the cc,us,ud families versus 
the experimental pairs of 3.526,3.525;1.402,1.375;1.231,1.303. The agreement of the theoret- 
ical and experimental mass differences is excellent for the ip system, slightly too large and 
of the wrong sign for the K system and too small and of the wrong sign for the ud system. 
Part of the cause of this pattern is that pure scalar confinement worsens the fit for the light 
mesons because of its tendency to reverse the spin-orbit splitting, thereby shifting the center 
of gravity. The agreement, however, for the hght systems is nevertheless considerably better 
than that in the case of the fine structure splitting R ratios. Another part of the discrepancy 
may be due to the uncertain status of the light ^Pq meson as well as the spin-mixing in the 
case of the K* multiplet. Note that in the case of unequal mass P states, our calculations 
of the two values incorporate the effects of the L ■ {si — S2) spin-mixing effects. (The use of 
nonrelativistic notation is only for convenience.) 

These differences between heavy and light meson systems also occur in the mixing due 
to the tensor term between radial S and orbital D excitations of the spin-triplet ground 
states. This mixing occurs most notably in the cc, us and ud systems. The three pairs of 
values that we obtain are 3.808,3.688;1.985,1.800;1.986,1.775 respectively versus the data 
3.770,3.686;1.714,1.412;1.700;1.450. Our results are quite reasonable for the charmonium 
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system but underestimate considerably the splitting for the light quark systems. As hap- 
pened for the significant disagreement in the case of the fine structure, our results here 
worsen significantly for the light meson systems. The spectroscopy of the lighter mesons is 
undoubtedly more complex due to their extreme instability (not accounted for in our ap- 
proach). Note, however, that for the spin-spin hyperfine splittings of the ground states the 
more relativistic (lighter quark) systems yield masses that agree at least as well with the 
experimental data as do the heavier systems. This same mixed behavior shows up again for 
the radial excitations. 

The incremental contributions for the six ^5"! states of the T system is just 1.8. It is 
12.9 over three states for the triplet charmonium system (primarily due to the ip deviation), 
3.0 for the two (p states, 1.6 for the three ^5*0 states of the K system (note, however that 
these fits include expected errors due to the lack of level shift mechanisms and are thus 
reduced), 7.3 for the two ^Si states of for the K* system, 2.2 for the three triplet ud states 
and 8.2 for the three singlet states . The contribution at first increases, then decreases 
with the hghter systems. Overall, the masses are much too large for the radially excited 
light quark mesons. These discrepancies may be due both to neglect of decay-induced level 
shifts and to the increased confining force for large r from linear to quadratic (there is no 
term to compensate for the quadratic S"^ term). 

The isospin splitting that we obtain for the spin singlet B meson system is 1 MeV. Our 
calculation includes the contribution from the u — d mass difference of 3.3 MeV as well as 
that due to different charge states. The effect of the latter tends to offset the effects of 
the former since the b and u have the same sign of the charge while the b and d have the 
opposite while > niu- In the experimental data this offset is complete (0). In the case of 
the D+ - splitting our mass difference of 7 MeV represents the combined effects of the 
u — d mass difference and the slightly increased electromagnetic binding present in the case 
of the and the slightly decreased binding in the case of the D'^. The experimental mass 
difference is just 4 MeV. These effects work in the same way for the spin-triplet sphtting 
resulting in the theoretical value 5 MeV compared with the experimental value 3 MeV. 
For the isodoublet we obtain 4 MeV versus about for the experimental values. Our 
isospin splittings are enhanced because of the large u — d quark mass difference that gives 
the best overall fit. For the K — K* family the experimental value for the isospin splitting 
is 4 MeV for the singlet and triplet ground states. This sphtting actually grows for the 



23 



orbital excitation {K2) to 7 MeV. The probable reason for this increase is that at the larger 
distances, the weak influence of the Coulomb differences becomes small while only the actual 
u — d mass difference influences the result (although it does seem rather large). It is difficult 
to understand why our results stay virtually zero for all three isodoublets. Note that as 
with the B doublets, the theoretical contributions of the combined effects of the u — d mass 
differences and the electrostatic effects tend to cancel. However, the experimental masses 
do not show this expected cancellation. 

1. Implications of our Model for the New 2.32 GeV D* Meson 

Recently, the BaBar Collaboration [58] found evidence for a new 0+ strange-charmed 
meson at 2.32 GcV. Using the parameters above and assuming the state is a ^Pq cs meson 
we find a predicted mass of 2.35 GeV, about 130 McV below our predicted value for the 
counterpart. The corresponding mass difference in the Godfrey-Isgur model is 2.590-2.480 
=110 MeV. Both are well off the experimental mark of 2.572-2.317=255 MeV. It is not 
surprising that its place in the quark model has been the subject of some debate. 

Overall comparison with the experimental data shows that the primary strength of our 
approach is that it provides very good estimates for the ground states for all families of 
mesons and for the radial excitation and fine structure splittings for the heavier mesons. 
On the other hand, it overestimates the radial and orbital excitations for the light mesons. 
Its worst results are those for the fine structure splittings for the us, ds and ud mesons. 
Both weaknesses are probably due to long distance scalar potential effects. Below, we shall 
discuss other aspects of our fit to the spectrum when we compare its results to those of other 
approaches to the relativistic two-body bound state problem. 

B. Explicit Numerical Construction of Meson Wave Functions 

There are 89 mesons in our fit to the meson spectrum. An important advantage of the 
constraint formalism is that its local wave equation provides us with a direct way to picture 
the wave functions. As examples, we present the wave functions that result from our overall 
spectral fit for three mesons: the tt (Figure 1) , for which we present the radial part of 
0+ = ■^i -I- ■04 that solves Eq.(A.5); 
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FIG. 1: The tt wave function plotted against x = log(r/ro) 




FIG. 2: The p{S) wave functions plotted against x = log(r/ro) 



the p (Figures 2 and 3) for which we present the radial parts of the wave functions 0+ 
for both S and D states that solve Eqs.(A.ll,A.17) ; and the ip/J (Figures 4 and 5) for 
which we present the radial parts of the wave functions 0+ for both S and D states that also 
solve Eq.(A.ll,A.17). In each plot the scale Tq is proportional to the Compton wavelength 
corresponding to the nonrelativistic reduced mass /i of the two quark system. In the table 
below, for each of the plotted mesons, we give the scale factor tq and the root mean square 
radius (in Fermis) computed from these meson wave functions. For the p and ip mesons we 
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FIG. 3: The p{D) wave function plotted against x = log(r/ro) 




also give the computed probabilities for residing in the S and D states. 

Meson tqH V< > S D 
TT 0.0004 0.21fm 1.00 0.0 
p 0.013 0.73fm 0.861 0.139 
ip 0.084 0.36fm 0.9974 0.0026 

Using a scheme outlined in Appendix B, we obtain an analytic approximation to the 
meson wave functions in terms of harmonic oscillator wave functions. The two primary 
parameters we use for each meson are the scale factor a and the leading power (short 
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FIG. 5: The V'(-C') wave function plotted against x = log(r/ro) 



distance behavior) exponent k. In addition we take as parameters the coefficients of the 
associated Laguerre polynomials. We write the radial wave function for each meson in the 
form 

N 

) = J^CnVn{r) (4.9) 



N 

u(r] 

n=0 

where 



\/ (n+f-l/2)! -"-^'/2)/LrV.(,.) 



""('■) ,: _j,r... <'M-yV2)y''Lt'/\y') (4.10) 

in which y — r/a — ae^ and (with z — y^) 

„z„-fe+l/2 Jn 

^ ' n\ dz^^ ' ^ ' 

We then vary the two parameters a and k to obtain the best fit. The coefficients are fixed 

by 

/• + 0O 

Cn= Vn{r)u{r)dr (4.12) 
Jo 

For meson radial wave functions with more than one component (like the i/j/ J) we fit each 
component separately. In the table below we give a typical list for parameters a, k, Cn for 
the TT, p, and ip/ J. 
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TT 


P 




k 


2.30734E-001 


9.85790E-001 


9.27248E-001 


0? 


1.22106E-004 


2.04708E-001 


5.85947E-002 


Co 


-9.70613E-001 


5.68290E-001 


8.63401E-001 


Cl 


1.97188E-001 


-5.54267E-001 


-3.77851E-001 


C2 


-1.18926E-001 


4.55647E-001 


2.70111E-001 


C3 


3.93232E-002 


-2.95969E-001 


-1.44888E-001 


C4 


-4.74935E-002 


2.11945E-001 


1.05621E-001 


C5 


1.59519E-002 


-1.29901E-001 


-5.85549E-002 


C6 


-2.21638E-002 


8.87707E-002 


4.46522E-002 


C7 


9.35388E-003 


-5.36537E-002 


-2.44101E-002 


C8 


-1.12997E-002 


3.57731E-002 


1.98781E-002 


Cg 


5.74799E-003 


-2.16185E-002 


-1.03167E-002 


ClO 


-6.24195E-003 


1.42167E-002 


9.24913E-003 


Cll 


3.44862E-003 


-8.57381E-003 


-4.34130E-003 


Cl2 


-3.63673E-003 


5.67698E-003 


4.49675E-003 


Cl3 


2.04307E-003 


-3.31349E-003 


-1.77086E-003 


Cl4 


-2.16019E-003 


2.33901E-003 


2.29266E-003 


Cl5 


1.22870E-003 


-1.19431E-003 


-6.63516E-004 


Cl6 


-1.26919E-003 


1.03806E-003 


1.23170E-003 


Cl7 


7.72030E-004 


-3.42741E-004 


-1.93158E-004 


Cl8 


-7.16255E-004 


5.20857E-004 


6.97788E-004 


Cig 


5.18700E-004 


-5.02603E-006 


-3.64677E-007 


C20 


-3.71156E-004 






C21 


3.77233E-004 






C22 


-1.56718E-004 







(4.13) 



We note several features. First, the fit to the tt wave function appears to converge signif- 
icantly more slowly than those for the p and i)jJ. (We do not present plots comparing 
the numerical wave function with the harmonic oscillator wave function fits since there are 
no visible differences). Also note that the tt's short distance behavior is distinctly different 
from those of the other two, having a stronger radial dependence at the origin. All three 
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wave functions possess polynomial coefficients that exhibit an oscillatory behavior. 
C. Numerical Evidence for Goldstone Boson Behavior 

In our equations, the pion is a Goldstone boson in the sense that its mass tends toward 
zero numerically in the limit in which the quark mass numerically goes toward zero. This 
may be seen in the accompanying plot Figure 6 (units are in MeV) . Note that the p meson 
mass approaches a finite value in the chiral limit. This non-Goldstone behavior also holds 
for the excited pion states. None of the alternative approaches discussed in the following 
sections have displayed this property. Another distinction we point out is that our u and 
d quark masses (on the order of 55-60 MeV) are significantly smaller than the constituent 
quark masses appearing in most all other models (on the order of 300 MeV) - closer to the 
small current quark masses of a few MeV. Note, however, that the shape of our pion curve 
is not what one would expect from the Goldberger-Trieman relation 

m^^mlF^. (4.14) 

Thus this aspect of our model requires further investigation. 

V. COMPARISON OF STRUCTURES OF TWO-BODY DIRAC EQUATIONS 
WITH THOSE OF ALTERNATIVE APPROACHES 

So far, we have obtained spectral results given by our equations when solved in their own 
most convenient form. In Sections (VI-IX) we shall compare our results with recent universal 
fits to the meson spectrum produced by a number of other authors. These approaches employ 
equations whose structures (at first sight) appear radically different from ours. However, 
as we have shown elsewhere [35], because our approach starts from a pair of coupled but 
compatible Dirac equations, these equations can be rearranged in a multitude of forms all 
possessing the same solutions. Among the rearrangements arc those with structures close 
to those of the authors whose spectral fits we shall shortly examine. In order to see how 
structural differences in each case may lead to differences in the resulting numerical spectra, 
we shall begin by considering relevant rearrangements of the two-body Dirac equations. 

The first two alternative approaches which we shall discuss use truncated versions of the 
Bethe-Salpeter equation (Salpeter and quasipotential) while the third uses a modified form 
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pi and rho masses vs. quark mass 




FIG. 6: TT and p masses versus quark mass in MeV 



of the Breit equation. In order to relate the detailed predictions of our approach to these 
alternatives, we need to relate our minimal substitution method for the introduction of in- 
teractions to the introduction of interaction through the use of kernels that dominates the 
older approaches. The field-theoretic kernel employs a direct product of gamma matrices 
times some function of the relative momentum or coordinate. What is the analog of the 
kernel in our approach? In earlier work we found that we could obtain our "external poten- 
tial" or "minimal interaction" form of our two-body Dirac equations from yet another form 
displaying a remarkable hyperbolic structure. We were able to recast our compatible Dirac 



30 



equations (2.1a,2.1b) as 

Sii; = (cosh(A)Si + sinh(A)S2)V' = 0, 

52V' = (cosh(A)S2 + sinh(A)Si)V' = 0, (5.1) 

in which [62] 

Slip = {Sio cosh(A) + sinh(A))^/' = 0, 

82^^ = (520 cosh(A) + cSio sinh(A))V^ = 0, (5.2) 
in terns of free Dirac operators 

Sio'ip = ( - Pi^i ■ P + ei/3i75i + mi-f5i)ijj 

= {(^2^2 ■ P + £2/^2752 + "^2752)^ (5.3) 

and the kernel 

A = ^751752 [-C'(a;±) + 7i • 72^(a;±)]- (5.4) 

We then recover the exphcit "external potential" forms of our equations, (2.1a,2.1b) from 
(5.1,5.2) by moving the free Dirac operators <Sjo to the right to operate on the wave function. 
This rearrangement produces the derivative recoil terms apparent in Eqs.(2.1a,2.1ba)). A 
may take any one of (or combination of) eight invariant forms. In terms of 

Oi = -751752, (5.5) 

these become A{x±) = -L{x±)Oi/2,-fi ■ P72 • Pj{x±)Oi/2,-fi± ■ -f2±G{x±)Oi/2 or ai ■ 
a2J-{x±)Oi/2 for scalar, time-like vector, space- like vector, or tensor (polar) interactions 
respectively. Note that in our A(x_l) in Eq.(5.4) above, Q{x±) enters multiplied by the 
electromagnetic-like combination 71 • 72 = —71 • -P72 • P + 7i± • 72± of time and space-like 
parts. This structure appears as a result of our use of the Lorentz gauge to introduce vector 
interactions in the classical version of the constraint equations or as a result of our use of 
the Feynman gauge to treat the field-theoretic version[55]. The axial counterparts to the 
constraints with polar interactions are given by (note the minus sign compared with the plus 
sign in Eqs.(5.1)) [62] 

Slip = (cosh(A)Si - sinh(A)S2)^/' = (5.6) 
Sii^ = (cosh(A)S2 - sinh(A)Si)V' = 0, 



31 



in which Si and S2 are still given by (5.2) with axial counterparts to the above A's given 



by C(xi)/2, 75171 • P75272 • PH{x±)Oi/2,-f5i-fi± ■ -f52-f2±I{x±)Oi/2 and ai ■ a2Y{x±)Oi/2 



respectively. The advantage of the hyperbolic form is that with its aid we may first choose 
among the 8 interaction types in an unambiguous way to introduce interaction (without 
struggling to restore compatibility) and then, for computational convenience, transform the 
Dirac equations to "external potential" form. In the weak-potential limit of our equations, 
the coefficients of 751752 in the expansion of our A interaction matrix in Eq.(5.4) directly 
correspond to the interaction kernels of the Bethe-Salpeter equation. Note however, that 
because of the hyperbolic structure, what we call a "vector interaction" actually corresponds 
to a particular combination of vector and pseudovector interactions in the older approaches 
(see Eq.(5.15) below). 

This difference in classification of interactions becomes apparent when we put our equa- 
tions into a Breit-like form. Consider the linear combination 



After simplification, the linear combination Eq.(5.7) of our two hyperbolic equations becomes 



/3i75iSi + /32752S2 



(5.7) 



For later convenience, form the interaction matrix 



Vix^) = -/5i75i/52752A(a;x). 



(5.8) 



= [Hio + H20 + V{x±, «!, a2, Pi, P2, 751, 752)]* 



(5.9) 



in which 



* = exp{-'D)'ijj 



(5.10) 



and 



Hio ^ ai-p± + Piuii, H20 = -a2 ■ p± + f32m2. 



(5.11) 



For the electromagnetic vector kernel A(x_l) = |[75i752]7i • 72^(a;±), T> then becomes 




(5.12) 



so that the relativistic Breit-like equation takes the cm. form 



= [ai • p - a2 • P + Pimi + P2m2 + w{l - exp[^(r)(ai • ot2 - 1)])]* 



(5.13) 
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In lowest order this equation takes on the famihar form for four- vector interactions (seem- 
ingly missing the traditional Darwin interaction piece 
~ f • aif • a2). 



= [ai • p - a2 • p + Ami + (32m2 - wQ{r){ai ■ a2 - 1)]*. 



(5.14) 



However, as we first showed in [45], expanding the simple structure of Eq.(5.13) to higher 
order in fact generates the correct Darwin dynamics. As a consequence, our unapproximated 
equation yields analytic and numerical agreement with the field theoretic spectrum through 
order a"^. Explicitly, our full interaction is 



so that our Breit-like potential contains a combination of "vector" and "pseudovector" in- 
teractions originating from the four- vector potentials of the original constraint equations in 
"external-potential" form. [63] 

In this section we have seen how the two-body Dirac equations with field-theoretic inter- 
action structure automatically retain the correct Darwin structure of QED. Such a demon- 
stration should be carried out for each alternative treatment (if possible) in order to check 
that truncations and numerical procedures have not destroyed its own version of the field- 
theoretic Darwin structure for its treatment of the vector interaction of QED (and associated 
vector structures in QCD). Explicitly in our own work we find that including all the cou- 
plings to smaller components of the wave function is crucial not only for our nonperturbative 
QED spectral results (see [25]) but also for our good results for tt — p splittings and the Gold- 
stone behavior of the pion as the quark mass tends toward zero. Without those couplings 
the good results for the positronium splittings and light mesons evaporate. 



exp[(ai ■ CX2 - l)Q] 



= [3cosh(^) + cosh(3^) + 751752(3 sinh(^) - sinh(3^)) 

+ ai • a2(sinh(3e?) + sinh(^)) + <ri • <r2(cosh(^) - sinh(3^))] (5.15) 
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VI. THE WISCONSIN MODEL OF GARA, DURAND, DURAND, AND NICK- 
ISCH 



A. Definition of The Model and Comparison of Structure with Two-Body Dirac 
Approach 

The authors of reference [48] base their analysis of quark- antiquark bound states on the 
reduced Salpeter equation containing a mixture of scalar and vector interactions between 
quarks of the same or different flavors. When rewritten in a notation that aids comparison 
with our approach, their bound state equation takes the cm. form 

[w-uj^-uj2Mp)=A^iph' J |^[^(p-p')7m^(p')7"+'5(p-p')$(p')]7°A-(-p) (6-1) 

in which A and S are functions that parametrize the electromagnetic-like and scalar inter- 
actions, A='= are projection operators, w is the cm. energy, coi = (p^ + m^)^/^, while $ is a 
four by four matrix wave function represented in block matrix form as 



(6.2) 



They obtain this equation from the full Bethe-Salpeter equation by making an assumption 
equivalent to using a position-space description in which they calculate the interaction po- 
tential with the equal time constraint, neglecting retardation. (These are the usual ad-hoc 
assumptions that in our approach are automatic consequences (in covariant form) of our 
two simultaneous, compatible Dirac equations.) These restrictions turn Eq.(6.1) into the 
standard Salpeter equation. In addition the Wisconsin group employs what we call the 
"weak potential assumption": ('u;-|-a;i-|-a;2) >> V. This assumption turns Eq.(6.1) into the 
reduced Salpeter equation which, because of the properties of the projection operator, allows 
the Wisconsin group to perform a Gordon reduction of its equation to obtain a reduced final 
equation in terms of 4>'^'^ alone. In our approach we make no such "weak potential assump- 
tion" and therefore must deal directly with the fact that our Dirac equations themselves 
relate components of the sixteen component wave function to each another. Unlike what 
happens in the reduced Salpeter equation, in our method this coupling leads to potential 
dependent denominators, a strong potential structure that we found crucial in demonstrat- 
ing that our formalism yields legitimate relativistic two-body equations. Just as we do, 
however, the Wisconsin group works in coordinate space where the dynamical potentials are 
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local and easy to handle. However, in their method upon Fourier transformation the kinetic 
factors cui then become nonlocal operators. In contrast, the entire dynamical structure of 
our two-body Dirac equations is local as long as the potentials are local. 

The Wisconsin group uses local static potentials that play the role of our Adler-Piran 
potential: 

Aa (r) , 3 
A{r) = --^^e-^' + S{-^ + Ar){l - e-^') 

(j V v 

S{r) = (1 - 5){-^ + Br){l - e'^') + (C + Cir + C2r')(l - e-'^'")e-'^" (6.3) 

Note that Gara et al introduce a confining electromagnetic-like vector potential propor- 
tional to a parameter 5. This differs from our approach in which the (dominant) linear 
portion of the confinement potential has no electromagnetic part. Like Adler's poten- 
tial, theirs has a long range 1/r part (the so-called Luscher term). Its short range part 
is electromagnetic-like just as is ours, and like Adler's is obtained from a renormalization 
group equation. 

They base their analysis on a nonperturbative, numerical solution of the reduced Salpeter 
equation Eq.(6.1) with interaction Eq.(6.3). 



B. Comparison of Wisconsin Fit with that of Two-Body Dirac Equations 

In Table II we include the Wisconsin variable-5 (vector and scalar confinement) best fit 
results, and the best fit results our method gives when restricted to the 25 mesons they 
consider. For uniformity of presentation we give all of the Wisconsin results in terms of 
absolute masses (rather than the mass differences and averages these authors presented for 
the spin-orbit triplets). Although Gara et al. did not perform the same fit that we do, 
we present (in parentheses) the incremental contribution for each meson so that we can 
easily compare the results of the two methods. We also compare their R values and avg. 
to ours directly in the discussion below. 

Our results are closer to the experimental results for 16 of the 25 mesons. In detail, 
their R values for the T and ip families of 0.83,0.78, and 0.60 are less accurate than two 
of our values of 0.64,0.68, 0.35 respectively. Their averages [5{^P2) + 3(^Pi) + l(^Po)]/9 
of 9.902 ,10.262, 3.513 and ours (9.901, 10.264, 3.513) are essentially the same compared to 
the experimental results of 9.900, 10.273, 3.525 MeV. Their hyperfine splittings for the two 
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charmonium multiplets of 200 and 47 MeV are significantly worse than our fits of 150 and 
79 MeV. Their hyperfine sphttings for the mesons with one d or s quark are 27, 51, and 127 
MeV. Our fits of 128,138, 420 MeV respectively are much closer to the experimental results 

of 141, 141, 398 MeV. 

The radial excitation energies for the two lowest T excitations and the singlet and triplet 
charmonium excitations are again accounted for significantly better by three of four of our 
values of 569,335,636,568 MeV for the results in the last column than by the Wisconsin 
results of 602,331,654,491 MeV. In summary, the major strength of our approach is reflected 
in its better fits to the hyperfine splittings and radial excitations. The Wisconsin group's 
results for the fine structure splitting are overall about the same as ours. Moreover, even a 
casual glance at the results shows one glaring discrepancy that results from their approach 
- their hyperfine splittings for the light quark mesons. The cause of this is probably the fact 
that their reduced Salpeter approach does not include coupling of the upper-upper piece to 
the other 12 components of the 16 component wave function. In fact, the lighter the meson, 
the worse is their result. In our QED numerical investigations we found that couplings to 
the other components of the wave function were essential in order to obtain agreement with 
the standard pcrturbativc spectral results of QED. Wc have found that the same strong- 
potential effects that led to our successful results in QED are responsible for the goodness 
of our hyperfine splitting, particularly for the mesons containing the light quarks. It would 
be important to test the Wisconsin group's procedure (with its deleted couphngs to the 
other wave functions) numerically with A — —a/r and 5" = for positronium to determine 
whether the problems that the Wisconsin model has with mesonic hyperfine splittings in 
QCD arc reflected in its results for QED. 

Gara et al. point out that in their approach the straight line Regge trajectories {j ver- 
sus w'^) for the light quark systems are much too steep, with slopes greater than twice the 
observed slopes for pure scalar confinement. The best fit experimental slope and intercept 
values for the p, a2,P3 trajectory are (0.88,0.48). The slope and intercept values that we 
obtain for our model in Table I are (0.87,0.47), in excellent agreement with the best exper- 
imental fit. For the 0i,/2,03 trajectory the experimental values arc (0.83,0.11) while our 
model of Table I produces the set of values (0.85,0.095). The intercepts are not as accurate 
as those for the p trajectory although our results actually produce a tighter fit to a straight 
line trajectory than do the experimental results. Finally we come to the tt, bi, 7r2 trajectory. 
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We obtain the values (0.57,-0.04). Compared to the experimental values of (0.72,-0.04) our 
slopes are about 25-30% small, although our fit to the straight line is just as tight. The 
probable reason for the relative advantage of our results over those of the Wisconsin group 
is that our bound state equations include a strong potential structure, and are not limited 
by the weak potential approximation built into the reduced Salpeter equation. 

VII. THE IOWA STATE MODEL OF SOMMERER, EL-HADY, SPENCE AND 
VARY 

A. Definition of The Model and Comparison of Structure with Two-Body Dirac 
Approach 

The Iowa State group introduces a new relativistic quasipotential reduction of the Bethe- 
Salpetcr equation. They use the well known fact that there are an infinite number of such 
reductions [65] to construct a formal quasipotential parametrized in terms of two indepen- 
dent constants. They show that when all of the most often used reductions are specialized to 
QED, they fail to numerically reproduce the correct ground state result for singlet positron- 
ium through order q;^[66]. These authors then fix the free parameters in their quasipotential 
by requiring that their resulting ground state energy lie close to the well-known perturbative 
value. In addition, the form of the quasipotential reduction they use produces a projection 
to positive energy states only. The Iowa State group uses a scalar linear confinement plus 
massless vector boson exchange-potential with the kernel 



The QCD coupling ag that they use is treated as a running coupling constant that depends 
on the momentum transfer and two parameters. Their quasipotential reduction incorporates 
zero relative energy in the cm. frame. 

B. Comparison of Fit with that of Constraint Approach 

In Table III, we give the Iowa State group's results for a set of mesons together with our 
results for the same set of mesons. In the fourth column of this table we present the results 
we would obtain from our approach if we limited our fit just to the 47 mesons used by the 




+ 47rb 



lim r d .2 7o X 7o 



(7.1) 



// — > '-9//-' —{q — q'y + fj,'^ 
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Iowa State group. We use the same RMS fitting procedure used by these authors instead of 
the fit used in our Table I. The results are quite similar, 50 for the Ohio State model 
and 53 for our model. 

Of the 47 mesons in their table, our fits are closer to the data in 25. Thus, accord- 
ing to this crude measure there is no significant difference between the results of the two 
approaches. [67] We proceed now with a detailed comparison. Their R values for the two bot- 
tomonium and one charmonium multiplets are 3.25,1.09,1.09. Our R values of 0.70,0.74,0.44 
are considerably closer to the experimental ratios of 0.66 (0.61), 0.56(0. 61), 0.47. (We make 
no comparison for the three light quark multiplets {ss, su, ud) since the Iowa State Group 
did not calculate the ^Pq states. ) We note, however, that for the pairs of su and ud 
their results for - ^Pi sphttings are substantially better than our results. In particular, 
unlike our results, theirs do not have an inversion of the splitting. Our poor results for 
these splittings are likely due to a larger infiuence of the scalar than the vector portion 
of the spin-orbit interaction. Comparing their ^P averages [5(^P2) + 3(^Pi) -|- l(^Po)]/9 of 
9.859,3.497,1.433,1.015 GeV for the lowest lying spin-orbit multiplets hsted in the table with 
our values of 9.902,3.516,1.470,1.386 and the experimental results of 9.900,3.527,1.503,1.303 
GeV we see that ours are closer in each case to the experimental results. Wc sec also that 
for charmonium, our average is nearly equal to our ^Pi level while the Iowa State results are 
75 MeV higher than their ^Pi level. For the ud system, our average is 25 MeV higher than 
our ^Pi level while theirs is 122 MeV above their calculated ^Pi level. Their values of the 
hyperfine {^Si Sq) splittings are 98, 48, 100,108,421,677 MeV for the two charmonium 
multiplets, and the D — D* , Dg — D* , K — K*, tt — p pairs. Comparison with the experimen- 
tal splittings of 117,92,142,139,398,628 MeV and our results of 159,82,137,156,376,593 MeV 
show the constraint results closer to the experimental splittings on all but the ground state 
charmonium pair. (We have commented earlier on the origin of the descrepancy between 
our value and the experimental result.) We next wish to compare the results generated 
in both approaches for the spin-spin effect embodied in the ^Pi - ^Pisphttings. For the 
cc, su, pairs the Iowa State results are 10,43,4 MeV compared to the experimental results 
of 15,136(129),28 (0) McV and the two-body Dirac results of 16,78,279 MeV. For the 
heavier two pairs, the constraint splitting results are substantially closer to the experimen- 
tal results. This resembles the similar spin-spin pattern found in the 5"— state hyperfine 
splittings. Our poor result for the ud meson has the same origin as our poor result for the 
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R value mentioned above. Finally, we compare the radial excitations. The six upsilon states 
in the experimental column of the table occur at intervals of 562,332,225,285,254 MeV. The 
three charmonium triplet states and the two charmonium singlet states occur at intervals of 
589,354,614 MeV while the two ss and ud states occur at intervals of 661 and 1160 MeV. 
The corresponding Iowa State intervals are 544,335,270,259,226,597,416,647,625,1304 MeV 
while our intervals are 578,345,260,218,191,560,395,637,753,1331 MeV. The Iowa State ra- 
dial excitation splittings are closer to the experimental values on four of the five upsilon 
splittings, one of the three charmonium splittings and both of the lighter quark splittings. 

Even though the RMS values obtained in each approach are nearly the same, on most 
of the detailed comparisons made above the constraint approach appears to give better 
fits. The exceptions to this arc the radial excitations and some of the heavier light-meson 
excitations. The largest portion of our RMS values come from the heavy-light meson orbital 
and radial excitations. 

We have long argued that any proposed relativistic wave equation should be tested in 
terms of its ability to reproduce known perturbative results of QED and other relevant 
relativistic quantum field theories when solved nonperturbatively before being applied to 
QCD. The Iowa State group in fact adopts this philosophy in order to resolve an ambiguity 
in the construction of the quasipotential in their wave equation by demanding that it re- 
produce the ground state level of singlet positronium numerically. This requirement fixes 
the values of the two parameters of their quasipotential mentioned above. In contrast, the 
constraint approach has no free parameters of the type used by [49] for the quasipotential 
reductions. Instead, its Green function is fixed. While within the constraint approach 
the connection between the kernel and the invariant constraint functions (e.g. Q, L) does 
involve some freedom of choice (see Eqs. (2. 14, 2. 15, 2. 17)), that freedom is not determined by 
the requirement that the model fit a particular state but instead is fixed by fundamental dy- 
namical requirements following equivalently from classical or from quantum field theory and 
resulting in the appearance of a minimal form of the potential (see Eq.(2.21) and below). 
Several features separate the two approaches. First, as we found in [25] the QED results 
provided by our equation agree with those of standard perturbative QED for more than just 
the ground state while it is unknown if the parameters that the Iowa State model uses that 
ensured its fit to the singlet ground state of positronium would work for the other states. 
Second the constraint approach generates similar structures for scalar interactions and sys- 
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terns of vector and scalar interactions with agreement with the corresponding perturbative 
field-theoretic results while again it is unknown whether the parameters that the Iowa State 
model uses that gave good fits to the singlet ground state of positronium would work in the 
presence of other potentials. Third, the match to singlet positronium that we obtained was 
an analytic consequence of our equations for QED and therefore a test of those equations 
[68], not the result of a numerical fit. Fourth, our approach includes essential contributions 
from all sixteen components of the relativistic wave function, not just the "positive energy" 
components [69]. Fifth, an important consequence of the fully relativistic dynamics and 
gauge-theoretic structure of the constraint equations is that they produce values of the light 
quark masses closer to current algebra values than do alternative approaches. The quark 
masses that we obtained in our comparison fit with the Iowa State model are rris = 314 
MeV and = — 67 MeV which are significantly closer to the current algebra values 
of 777.5 ~ 125 MeV and 777^, tti^ ~ 3 — 6 MeV than the Iowa State model's values of 405 and 
346 MeV respectively. 

YIII. THE BREIT EQUATION MODEL OF BRAYSHAW 

A. Definition of The Model and Comparison of Structure with Two-Body Dirac 
Approach 

Brayshaw [50] treats quarkonium with the aid of the Breit equation and an interaction 
Hamiltonian with five distinct parts, four of which are independent. As usually done for 
the Breit equation the times associated with each particle are identified or related in some 
favored frame (normally cm.) selected so that the relative time does not enter the potential. 
In that frame Brayshaw uses the equation 



= {Ho + Hc + Hb + Hs + Hi + Hl)-^ = 



(8.1) 



in which Hq is the free Breit Hamiltonian 
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while He and Hb are a Coulomb and an associated Breit interaction 




ci(ai • a2 + ai • ra2 ■ r) 

2t 



(8.3) 



As indicated in our discussion about the Salpeter equation in Section(VI), this part of the 
interaction comes from the vector portion of the kernel. The author acknowledges the dif- 
ficulties associated with the Breit interaction, pointing out that the radial equation has a 
singularity at a radial separation of tq = —c\jw > 0. He bypasses Breit 's proposal that 
this interaction be used only in first order perturbation theory by using only positive en- 
ergy spinors in his variational procedures. We point out that this was not necessary in 
our approach since the hyperbolic structure of our eight basic interactions avoids problems 
inherent in Breit 's formulation [9]. In particular, it avoids appearance of midpoint singular- 
ities. Unfortunately, just like the Wisconsin group, having avoided the pitfalls of the Breit 
equation, he uses his replacement without testing whether or not his formalism would yield 
the standard QED results numerically if he limited his interaction to the usual Coulomb 
interaction. Once again such a test would (if successful) help eliminate the possibility that 
the wave equation introduces spurious physics. 

In Eq.(8.1), Hl is a long range confining portion which incorporates the requirement that 
the wave function vanish identically for radial separations r > a with a boundary condition 
at r = a. Brayshaw argues for this term over and above a linear confinement piece on the 
grounds that at some separation Vp corresponding to a threshold energy Ep , production 
of qq pairs should become energetically favorable. His radial parameter a plays the role of 
Tp in specifying the range at which such effects (among others) dominate confinement. He 
expects that a is on the order of (r) for the light quark mesons while wave functions for the 
heavy quark mesons would have fallen to zero for r << a. When introducing the explicit 
form of his linear confinement potential, the author finds that it cannot simply be added as 
a Lorentz scalar to the Hamiltonian since such a term produces far too large a mass shift 
for the light quark systems. Instead he chooses 



which he shows contributes very weakly for the light quark systems, while contributing 
significantly for the heavy quark systems with an intermediate contribution for the hydrogen- 



Hi = C2{Pi + P2)r. 
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like intermediate mass mesons. Unfortunately, however, we note the important fact that the 
Lorentz transformation character of this confining interaction is ambiguous, being neither 
scalar (~ /3i/32) nor (time-like) vector (~ lil2)- 

Finally Brayshaw introduces a special short range attractive piece solely in order to obtain 
a good fit to the pion and kaon. Instead of a spin-dependent contact term used in a number 
of semirelativistic approaches [51, 70, 71] he uses 

Hs = Hsihh + W ^/''"^^^"'^^ , (8.5) 

2(mi + m2 + C4) 

This term resembles a cross term between a linear confinement piece and the Breit term 
that might emerge from some sort of iteration. The short range character of this part- 
scalar, part- vector interaction is specified through taking b « a. In contrast, our approach 
possesses a short range spin-spin interaction that is quantum mechanically well defined and 
which arises straightforwardly from the Schrodinger reduction of our Dirac equations. We 
do not need to add it in by hand. 



B. Comparison of Fit with that of Constraint Approach 

In spite of its ad hoc nature, we have included the procedure of Brayshaw among our 
comparisons because it turns out that his resultant fit for the 56 mesons (that overlap with 
our fit) is quite good, just shghtly worse than our fit. In Table IV we include in the fourth 
column the fit we would obtain with our model if we included only the 56 mesons that 
our fit has in common with Brayshaw's. On a meson by meson basis we compare by using 
incremental values. 

Of the 56 mesons in the table, our fits are closer to data in only 26, although overall our fit 
is better. However, this overall difference may not be as significant as in the previous exam- 
ples because here we did not use identical fitting procedures for both models. Brayshaw's it! 
values for the two upsilon, the one charmonium, the K*, and p — tt triplet P multiplets are 
0.47,0.34,0.32,0.55,0.25,0.19 and are distinctly different from our values of 0.66,0.69,0.39,- 
0.71,-0.25,-5.67 and the experimental numbers of 0.66,0.61,0.48,0.09,-0.97,-0.4. Although the 
constraint/ Adler-Piran combination is distinctly better than the Breit/Brayshaw approach 
for the heavier mesons, both give poor R results for the lighter mesons. All of his light spin- 
orbit multiplets have masses that increase monotonically with j, unlike the pattern of the 
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experimental numbers. Although our results show a non-monotonic pattern that pattern 
also differs from that of the data. Note that the details of our patterns are greatly influenced 
by the presence of the scalar potential. Brayshaw's approach includes (see Hs) a partial 
Hamiltonian that governs intermediate range behavior, in which time-like and scalar inter- 
actions contribute equally. This may be responsible for the difference between his montonic 
pattern and that displayed by the data. 

Comparing his averages [5(^P2) + 3(^Pi) + l(^Po)]/9 to the ^Pi mesons 

for the charmonium, K*, and p — tt systems we find the following three 
pairs of numbers: 3. 517,3.498;!. 335,1.355;!. 251, 1.202. Comparison to our num- 
bers of 3.519,3.520;1.435,1.421;1.434,1.411 and the experimental numbers of 
3. 526,3.525;!. 402,1.375;!. 231, 1.303 shows that our approach gives better agreement 
for the heavier mesons, his somewhat better for the lighter while both do about the same 
for the K*. 

His values of the hyperfine splittings are 118,100,143,158,410,636 MeV for the two char- 
monium multiplets, and the D — D*, Dg — D^, K — K*,7r — p pairs. Comparing with the 
experimental splittings of 117,92,142,144,398,627 MeV shows a clear pattern of excellent to 
good results for the heaviest, lightest, and the intermediate more hydrogen-like mesons. Our 
results are 151,79,133,145,416,647 MeV. Our ground state charmonium result is not nearly 
as good as Brayshaw's while for the others we have about the same quality of fit. It may be 
that his choice of Hs rectifies the problem our treatment encounters. But, the disadvantage 
of this is that his R values for the heavy mesons are worse. This effect appears to be similar 
to the trouble we encountered, mentioned in our discussion of Table I in Sec. IVA. 

For the radial excitations, the four upsilon states in the data portion of the table occur 
at intervals of 563,332,225 MeV while the three charmonium triplet states and the two 
charmonium singlet states occur at intervals of 589,354,614 MeV. The pion excitation is 1160 
MeV. The corresponding Brayshaw intervals are 555,335,320,551,566,569,888 MeV while our 
intervals are 572,337,257,564,395,636,1403 MeV. With the exception of the second radial 
triplet upsilonium and charmonium excitation intervals, the fits of both models are of about 
the same quality. Note that excited pion predictions bracket the experimental results. This 
appears to be a common feature of the radial and orbital excitations of the light quark 
mesons, with his results on average closer to the experimental values. Our results are, on 
average, better for the heavier mesons. 
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However, his apparently good fit emerges from a potential structure that has ambiguous 
Lorentz transformation properties. The potentials are chosen in a patchwork manner us- 
ing the 5 parameters a, Ci, C2, C3, C4(he sets b = a/10). In terms of Lorentz transformation 
properties his scheme uses four invariant functions (scalar, time-like, electromagnetic like 
and mixed {Hs,Hb and Hj)). The Adler-Piran potential that wc use has only two invari- 
ant functions corresponding to scalar and electromagnetic like interactions. The constraint 
approach is not a patchwork; instead its wave equation itself (once A and S are chosen) 
fixes the spin, orbital and radial aspects of its potential and its spectra. We also note that 
just as in the case of the Wisconsin model, Brayshaw has not tested the nonperturbative 
reliability of his equation. On the other hand an important result of his approach is that 
the u, d quark masses required for his fit are very small (10 MeV) and significantly closer to 
the current quark mass values than ours. His strange quark mass (200 MeV) is also closer 
to the proposed current quark mass values than our value. 

The most important warning provided by Brayshaw's approach is that an ad hoc structure 
with ambiguous Lorentz properties can do so well at fitting the spectrum. 

IX. THE SEMIRELATIVISTIC MODEL OF GODFREY AND ISGUR 

A. Definition of The Model and Comparison of Structure with Two-Body Dirac 
Approach 

We begin with a general discussion of Semirelativistic Quark Models (with and without 
full Relativistic Kinematics). We term a "semirelativistic quark model" one that uses a 
two-body wave equation that takes one of the following three forms in the cm. frame: 

(p^ $(r, 81,82))^' ^{w - mi- 1712)1^ 

(y^p2"+mf + \/p2 + 7712 + $(r. Si, S2))V^ = 

(p2 + $(r,8i,82))V' = 6'(w)V'. (9.1) 

In each of these equations p^ is the square of the cm. relative momentum while $(r, Si, 82) 
is an effective potential which includes central, spin-orbit, spin-spin, tensor and possibly 
Darwin terms. In each, the wave function has four components with no coupling to lower- 
lower components. The most important difference between the first form and the others is 



44 



that the latter two have exact relativistic kinematics. The former is almost always called a 
nonrelativistic quark model although strictly speaking almost all spin dependences (at least 
those that arise from vector and scalar interactions) vanish in the nonrelativistic limit. These 
equations differ from the Two-Body Dirac equations and the Breit and instantaneous Bethe- 
Salpeter approaches primarily in that their spin-dependences are put in by hand, abstracted 
from the Fermi-Breit reductions of the Breit and instantaneous Bethe-Salpeter approaches. 
For Coulomb-like potentials originating in the Coulomb Gauge, these terms contain singular 
potentials. Consequently they must either be treated purely perturbatively (thus ruling out 
application to the light quark mesons) or through the introduction of smoothing parameters 
that may or may not be features of the actual potential. The two-body Dirac equations 
of constraint dynamics, like their one-body cousin, have a natural smoothing mechanism - 
potential dependent denominators in the spin-dependent and Darwin terms of the resultant 
Schrodinger-like form - that eliminates the necessity for ad hoc introduction of such terms. 
The Breit equation may also possess a natural smoothing mechanism, but a nonperturbative 
treatment of it leads to erroneous results in QED [5] . The instantaneous Salpeter equation 
may have a natural smoothing mechanism, but has not been tested nonperturbatively for 
QED even though the equation is over 50 years old. Authors who have attempted to use these 
types of semi-relativistic equations to treat the entire meson spectrum include Lichtenberg 
[70] (the third type), Stanley and Robson [71] and Godfrey and Isgur [51] (the second type), 
and Morpurgo, Ono, and Schoberl [46] (the first type) . Each of these authors ignore the 
spin-independent part of the Fermi-Breit interaction. This neglect is not justifiable since 
this part of the interaction will have an effect on S states that is significantly different 
from its effect on non ^-states, being normally short ranged compared with the rest of the 
central force part of the problem. In this paper, we select one of these models for our final 
comparison, the model of Godfrey and Isgur, since this model, even though already 18 years 
old, is by far the most often cited in recent experimental works and theoretical papers on 
rival approaches. 

As we have said, Godfrey and Isgur assume a semi-relativistic wave equation of the second 
type possessing exact relativistic kinematics but through the inconvenient sum-of-square- 
roots form. They then determine the form of interaction in the following way. They assume 
that the confining piece of the interaction is a world scalar. They modify the Coulomb 
potential with the aid of a smoothing function. At the same time they appear to ignore 
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the Darwin term (e.g. the spin independent contact term present in the one-body hmit) 
in the on-shell reduction of the qq scattering amphtude. Although they modify the short 
range part of their interaction with the aid of a smearing function, this modification does not 
compensate for the ignored Darwin term. We have shown elsewhere [28] , [30] that the Darwin 
interactions for scalar and vector interactions lead, through a canonical transformation to 
the quadratic local terms S"^ and J\? that appear in our equations. Since the authors have 
ignored this part of the Darwin interactions their results contain none of the dynamical 
consequences of the S"^ or — J? pieces. What portion of the Darwin term they include they 
parametrize separately just as they do the other portions of the Fermi-Breit interaction. 
These terms include the spin-spin contact term, the spin-orbit terms, and the tensor terms. 
In our opinion, this patchwork way of handling the physics blurs the relativistic significance 
of their quark model. In our two-body Dirac equations the Darwin portion and each of 
the spin-dependent portions is tied directly to and fixed by the Lorentz forms L{x_i), G{x±) 
of the interaction which are in turn set by the A invariant potentials. In QED these 
fixed terms yield the correct spectrum with no additional parameters needed to adjust their 
relative sizes. 

In addition to bypassing the problems of singular spin-dependent terms by assuming a 
smoothing parameter, Godfrey and Isgur include nonlocal (momentum-dependent) poten- 
tials by replacing the mass dependent mj^ in the Fermi-Breit term by (p^ -|-m^)~^/^. They 
claim that this is necessary because the Fermi-Breit reduction (or the on-shell qq scattering 
amplitude in cm.) does not adequately express the full momentum dependence (or nonlo- 
cal nature) of the potential. While this might be true, we have found that such nonlocal 
behavior is not necessary to obtain very good results either in lowest order QED or in the 
quark model. 

Like the Adler-Piran potential that we use in our approach, their potential includes a 
running coupling constant. In fact, by convolving a parametric Gaussian fit to the running 
coupling constant with the ^ , they obtain their desired smoothing of the Coulomb potential, 
thus killing two birds with one stone. In addition, they are able to treat the zero isospin 
mesons like the 7] and rj' by including a phenomenological annihilation term. We leave out 
this term in our results of Table 1-IV and in our comparison with the results of Godfrey and 
Isgur in Table V. Lichtenberg [70] has compared an earlier version of our quark model for 
the meson spectrum with that of Godfrey and Isgur. The potential we used in that earlier 
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version was the one-parameter Richardson potential, with the confinement piece chosen to 
be one-half time-hke vector and one-half scalar. As Lichtenberg pointed out, Godfrey and 
Isgur obtained significantly better agreement with the data than we did. He states that this 
is because they use significantly more parameters than wc do including four in the potential 
and six to describe relativistic effects, ten altogether, compared to our one. However, we do 
not believe that as a general rule the number of parameters that appear in the potential is, in 
itself, of as much significance as how these parameters are distributed. For example, in our 
present and previous models there are two invariant functions, A and S related to the single 
nonrelativistic (Adler-Piran) Vap that itself depends on two parameters. These parametric 
functions are not entirely independent, being related by Eqs. (3. 1,3. 8a, 3. 8b). Specifying their 
form fixes both spin- independent and spin-dependent parts of the quasipotential We 
might say that our formalism has 5 quark mass parameters and two parametric functions. 
Increasing the number of parameters that A, S depend on may or may not increase the 
goodness of the fit. According to our way of counting, Godfrey and Isgur have independent 
parametric functions for the two spin-orbit parts of the potential, the spin-spin contact part, 
the tensor part, the scalar potential, and the spin- independent part of the vector potential, 
altogether 6 parametric functions. From our way of counting the number of parameters 
the number of parametric functions would not increase no matter how many parameters 
are included in fixing the functional form of each of these six functions. Likewise, in our 
case, no matter how many parameters we use in fixing A, S there are only two independent 
parametric functions. Our approach is distinct from that of Godfrey and Isgur in that we 
do not alter the functional form at the level of the spin-dependence but rather at the level 
of the kernels. 

Finally, before we compare our present work with that of Godfrey and Isgur, we note that 
our present model differs from our earlier one used by Lichtenberg in his comparison of the 
two approaches. Our present treatment differs in its replacement of the Richardson potential 
by the Adler-Piran potential. The intermediate range form of the A-P potential is closely 
tied to an effective field theory related to QCD and is therefore superior to Richardson's 
ansatz. Furthermore, in calculations based on our earlier treatment we ignored the tensor 
coupling and unequal mass spin-orbit difference couplings which we explicitly include in the 
present calculations. We have also corrected a defect in the decoupling we used between the 
upper-upper and lower-lower components of the wave functions for spin-triplet states in our 
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older treatment. 

B. Comparison of Fit with that of Constraint Approach 

We now compare the fit given by our model to that provided by the model of Godfrey 

and Isgur. 

In Table V we display in the fourth column the fit we would obtain with our model if 
we included only the 77 mesons that our fit has in common with that of Godfrey and Isgur. 
We then compare the fits by examining the incremental values for each meson. (In an 
RMS comparison they would obtain about 63 compared to our value of 79). 

For the 77 mesons in their table, our fits are closer to data in only 32; overall their fit 
is better. Generally speaking our results are better on the newer mesons while their fit is 
better on the older mesons. A detailed comparison reveals the following. Their R values 
for the two upsilon, the one charmonium, the K* ^ (f) and p — tt triplet P multiplets are 
0.29,0.50,0.57,0.36,0.42,0.47 and are distinctly different from our values of 0.68,0.76,0.41,- 
0.66,-0.21,-4.00 and the experimental numbers of 0.66,0.61,0.48,0.09,-0.97,-0.4. As was true 
for the Brayshaw analysis, the constraint/Adler-Piran combination gives a distinctly better 
fit than the Isgur- Wise approach for the heavier mesons, while both give poor results for 
the lighter mesons. As was the case for Brayshaw's spectrum, none of their light multiplets 
are inverted, whereas although ours are inverted they are not inverted in the same way as 
the experimental numbers are. Again, our inversions are due to the action of the scalar 
potential. Godfrey and Isgur include a time-like contribution in the spin-orbit part of their 
Hamiltonian. This may be responsible for their lack of the partial inversion that appears in 
the data. 

Computing their averages [5(^P2) + 3(^Pi) -|- l(^Po)]/9 along with the 

^Pi mesons for the charmonium, K* and p — n system we find the follow- 
ing three pairs of numbers: 3.524,3.520;1.392,1.340;1.262,1.220. Comparison with 
our numbers of 3.519,3.520;1.424,1.411;1.419,1.397 and the experimental numbers of 
3.526,3.525;1.402,1.375;1.231,1.303 shows the constraint approach giving slightly better num- 
bers for the heavier mesons and the K* while the Godfrey-lsgur results are somewhat better 
for the lighter mesons. Their average [7{^D2) + 5{^Di) + 3(^Di)]/15 and their 
meson for the K* are 1.795,1.780 MeV while our results and the experimental results are 
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1.873,1.879 and 1.774,1.773 MeV respectively. Our results are relatively closer to one an- 
other while theirs are closer to the data in an absolute sense. This is indicative of the general 
trend of our orbitally excited hght mesons being somewhat high. We suspect that this is 
due to the S"^ behavior becoming dominant at longer distance, changing the behavior of the 
confining potential in the effective Schrodinger-like equation from linear to quadratic. 

Their values of the hyperfine splittings are 130,60,160,150,430,130,620,150,120 MeV for 
the two charmonium multiplets, and the D — D* ,Ds — D*, two K — K*, and three tt — p 
pairs. Comparison with the experimental splittings of 117,92,142,144,398,-48,627,165,354 
MeV and our results of 150,78,133,145,403,208,645,239,166 MeV demonstrates that while 
our results are closer than theirs for most of the newer mesons and the K — K*, their results 
are more in line for most of the older mesons. Again this shows a pattern of our method 
overestimating the radially excited states of the light mesons. 

Let us see if this trend of overestimation by the constraint approach continues for 
the radial excitations of fixed quantum numbers. The six upsilon states in the data 
portion of the table occur at intervals of 563,332,225,285,154 MeV while the three 
charmonium triplet states and the two charmonium singlet states occur at intervals 
of 589,354,614 MeV whereas the three singlet K and the two triplet K* states oc- 
cur at intervals of 977,370 and 520 MeV. Finally the three pion and three rho excita- 
tions occur at 1160,495 and 698,654 MeV. The corresponding Isgur-Wise intervals are 
540,350,280,250,220,580,420,650,980, 570,680,1150,580,680,550 MeV compared to our inter- 
vals of 570,336,256,213,186,561,393,633,1099,495,894,1383,634,986,561 MeV. Again we en- 
counter a pattern of our results being more accurate overall for the newer mesons while 
theirs are more accurate for the older ones (with our results too large for all of the older 
ones) . 

Primarily what we learn from this comparison is that not only does the scalar interaction 
lead to partial triplet inversions for the lighter mesons but also yields radial and orbital 
excitations that are too high for a related reason: the presence of the S'^ term in the effective 
potential. On the other hand, as Godfrey and Isgur themselves point out, their treatment of 
the relativistic effects is schematic, with no wave equation involved, allowing an uncontrolled 
approach in which there are no tightly fixed connections among the various spin-dependent 
and spin- independent parts of the effective potential 

An important feature of our approach that differs significantly from the model of Godfrey 
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and Isgur (as well as those of the Wisconsin and Iowa State groups ) is the size of its resulting 
light quark masses. Our u, d quark masses are about a factor of four or five smaller than 
theirs, significantly closer to the current algebra values. Godfrey and Isgur argue that 
since a constituent quark model requires dressed quarks of a finite size (to avoid singular 
potentials in their wave equation among other reasons) one should not expect the model 
quarks to have current-quark masses. We argue that a properly structured relativistic wave 
equation should not require finite quark sizes. Similar remarks have been made historically to 
justify tampering with the wave equation in QED to avoid treating singular terms. However, 
in QED those terms are perturbative artifacts. In fact, in the constraint equations for QED, 
they arise from premature weak-potential approximation to terms that arc actually well- 
behaved at the origin. Similarly, when we apply the constraint approach to QCD we need 
no size parameters. 

Finally we mention what we consider the major theoretical shortcoming in the approach 
of Godfrey and Isgur. The formalism that they use gives very good results on the hyperfine 
splittings for light and heavy mesons. However, it is unknown if this is an artifact of their 
smearing factors and the introduction of relativistic momentum dependent corrections to 
the potentials (that is, through the replacement of quark masses m by \/p'^ + rn?) needed 
to modify the singular nature of the potentials that they start with. 

It would be of interest to test the wave equation used by Godfrey and Isgur numerically 
with A — —a/r and S — Q ior positronium to see if any of their successes with mesonic 
hyperfine splittings are refiections of corresponding nonperturbative successes in QED. If 
their method were not able to obtain an acceptable fit to the QED spectral results through 
order a^, then the legitimacy of its fits in QCD would be seriously called into question. 
Without such tests one could not be sure whether the method they employ to avoid the 
singular potentials has distorted the dynamics. The constraint approach has passed this test 
in that without introducing additional parameters it does faithfully reproduce the correct 
spectral results in QED. 
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X. CONCLUSION AND WARNINGS ABOUT THE DANGERS OF "RELATIVIS- 
TIC" AND "NONRELATIVISTIC" SPECTRAL FITS 



In this paper, we have investigated how well the relativistic constraint approach performs 
in comparison with selected alternatives when used to produce a single fit of experimental 
results over the whole meson spectrum. This approach is distinguished from others by its 
foundation - a set of coupled, compatible, fully covariant wave equations whose nonpertur- 
bative numerical solution yields the mass spectrum along with wave functions for the qq 
meson bound states. Its virtue - generation of fully covariant spin structures - also serves 
to restrict and relate plausible interaction terms just as the ordinary single-particle Dirac 
equation determines relations among Pauli spin dependences and fixes the proper strength of 
the Thomas precession term in electrodynamics. The dynamical structures of the constraint 
approach were originally discovered in classical relativistic mechanics but have since been 
verified for electrodynamics through diagrammatic summation in quantum field theory in 
the field-theoretic eikonal approximation [27]. 

To use such relativistic equations to treat the phenomenological chromodynamic qq 
bound-state, one must construct a relativistic interaction that possesses the limiting be- 
haviors of QCD. In our approach we have done this by using the nonrelativistic static 
Adler-Piran potential to construct a plausible relativistic interaction that regenerates the 
AP potential as its nonrelativistic limit. In our equations, this process generates a host of ac- 
companying interaction terms. When describing these interactions, one must guard against 
a semantic difficulty in the verbal classification of the various parts of the interaction as 
"scalar", "vector", "pseudovector" etc. The various formahsms classify these in different 
ways but in our equations, the meaning of these terms can be readily determined through 
examining their roles in the defining equation Eq.(2. 26, 5. 1,5.2,5.4). Once these terms have 
been introduced, the constraint formalism automatically produces a system of important 
accompanying terms like quadratic terms that dominate at long distance (reinforcing or 
undermining confinement) or spin dependences that accompany chosen interactions produc- 
ing level splits that agree or disagree with the experimental results in various parts of the 
spectrum. 

After identification of the relativistic transformation properties of interaction terms the 
constraint method leaves almost no leeway for fiddling with (unnecessary) cutoffs, etc. Some 
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years ago, when applied to the e~e^ system, its structure proved restrictive enough to rule 
out within it the presence of postulated anomalous resonances [2 5, 72]. In recent work on 
the relation of our equations to the Breit and earlier Eddington-Gaunt equations for electro- 
magnetic bound-states, the method has explicitly demonstrated the importance of keeping 
spin couplings among pieces of the full 16-componcnt wave-functions whose counterparts are 
often truncated or discarded in alternative treatments [9, 37]. 

The fits that we have examined as alternatives fall into different classes: motivated rel- 
ativistic fits ( constraint vs truncations of standard field-theoretic), ad- hoc relativistic fits, 
and cautious semirelativistic fits. 

Among the relativistic ones, there is a danger exemplified by the Brayshaw model which 
achieves relative success despite the dubious relativistic nature of its interaction. As always, 
what makes fits hard to judge at this stage is the ease with which one can achieve apparent 
success over limited regions of the spectrum using highly-parametrized interactions. We 
have attempted to avoid this problem by limiting comparisons to published treatments that 
include both the light and heavy meson portions of the spectrum, not just one of the two 
sectors. Our choices for comparison are meant to be representative (we do not attempt an 
exhaustive review) (sec [73] for other important treatments). 

With the exception of the Iowa State model [49] all of the comparison models fail to 
test whether or not a nonperturbative treatment of their wave equations would yield the 
known results if the QCD kernels used were to be replaced by ones appropriate for QED. 
With the exception of the quark masses obtained by Brayshaw, [50] our light quark masses 
are substantially closer to the current algebra values than are those produced by the other 
comparison models. In our application of the constraint approach, it is possible to describe 
the nonperturbative physics that accommodates a typical size for an effective or constituent 
quark mass used in the other approaches and which at the same time has the size necessary 
to account for baryon magnetic moments. Even though our u and d quark masses are small 
compared with constituent quark masses found in the competing approach, if we compute 
the expectation value (Mj(^, S)) we find a range that includes those values. We find the 
range of values for this effective mass from 64 MeV for the pion to 390 MeV for the rho. 
Its value depends not only on the quantum numbers of the meson but also the flavor of the 
other quark. For example, for the D meson we find {M^^A, S)) — 190 MeV whereas for the 
B we obtain 258 MeV. 



52 



Finally, some authors have even produced unabashedly nonrelativistic fits. They claim 
to obtain good fits to the meson spectrum through the use of variants of the nonrelativistic 
quark model (NRQM). [46], [74]. These authors even claim success at fitting the light 

quark mesons for which the assumptions T « mc^, \V\ << mc^ of the nonrelativistic 
Schrodinger equation are patently false. What can account for the apparent success of the 
NRQM? 

Morpurgo states [46] that the various potential models, including the nonrelativistic quark 
model, are merely different parametrizations of an underlying exact QCD Lagrangian de- 
scription. That is, all use essentially the same spin and fiavor structures. For example, 

for the mesons one can derive a "parametrized mass" with general form (for the present 
discussion restricted to tt, K, p, K*) 

"parametrized mass" = A + B{P^ + P|) + Ccr^ -a^ +D{P^ + P^)ai-a2 (10.1) 

in which Pf is the projector onto the strange quark sector. These authors say that this 
structure although typical of an NRQM description, follows from QCD itself. They state 
that the form Eq.(lO.l) is common to all of the relativistic or semirelativistic quark models. 
They assert that any one of them can be successful but not superior to any other, if it merely 
reproduces the spin fiavor structure of the general parametrization. Thus, from their point 
of view selection of the "best" model is entirely a matter of taste and simplicity. We disagree 
with this assessment for the following reasons. First, the kinetic and potential parameters 
have significances beyond simply producing a fit for the two-body bound-state sector in 
isolation. When the spin-fiavor structure in (10.1) appears in the constraint approach, its 
accompanying constituent quark masses turn out to be closer to the current-quark masses 
than those produced by most other approaches while the constraint method requires only 
two parametric functions to be used beyond the parameters of the constituent quark masses. 
The constraint scheme successfully uses one set of these parametric functions for the entire 
spectrum of meson states including the radial as well as orbital excitations. But most im- 
portantly, within the bound-state spectrum itself, in our relativistic approach even though 
superficially sharing the basic spin-fiavor structure (10.1), all potentials do not fare equally 
well. The essential point is that even in the simplest form of our equations, the parametriza- 
tion is different from that given in the Morpurgo form in that its parameters A,B,C,D, are 
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themselves dependent on the energy operator on the left hand side. When that happens, 
some relativistic potentials do better than others. In particular, of those we investigated, 
the potential that works the best (the Adler-Piran potential) is one possessing many of the 
features important in lattice QCD calculations (e.g. linear and subdominant logarithmic 
confining pieces). The combination of the constraint approach with the Adler-Piran poten- 
tial embodies more of the important physical effects contained in QCD-related effective or 
numerical field theories. 

Can one understand the apparent successes of the NRQM fits by starting from the rela- 
tivistic treatments? Some authors [75] and [76] have used bounds on the kinetic square-root 
operator a/p^ + to attempt to understand the apparent success of the iiom^elativistic 
potential models for relativistic quark- ant iquark states. Instead, we will give an explanation 
that starts directly from the relativistic constraint approach. 

Some years ago, Caswell and Lepage [77] rewrote a relativistic constraint equation in 
an effective nonrelativistic Schrodinger-like form. Here, we do the opposite and recast the 
NRQM Schrodinger equation in a form resembling the constraint equation. As we have seen 
our two-body Dirac equations lead to an effective Schrodinger-like equation of the form 

\p^ + ^U^±,ai,a2M^b\w)i^ (10.2) 

In the cm. system this becomes 

[p' + $^(r,(ri,o-2)]V = &'(«^)V' (10.3) 
Even though the stationary state nonrelativistic Schrodinger equation 

[^ + Vir,a,,(T2M = EB^ (10.4) 

has a similar form, the corresponding structures in each have entirely different physical 
significances. For example, in Eq.(10.4), the vectors p and r are nonrelativistic quantities 
in contrast with their counterparts in the constraint approach that appear in the relativistic 
equation in the cm. system. One can easily manipulate the nonrelativistic Schrodinger 
equation into a form similar in appearance to the constraint Schrodinger form by multiplying 
both sides of the equation by 2// and adding ^^{w) — 2ij,Eb to both sides. The result is 

[p2 + $^(r, o-i, (T2)]iJ = b\w)ij (10.5) 
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in which 



$^(r, <Ti, 0-2) = 2//y(r, o-i, (T2) + b'^iw) - 2hEb 



(10.6) 



In numerical calculations the p operator and r variable are treated in the same manner 
in calculations based on both the relativistic constraint equation and the nonrelativistic 
equation. But as we have seen, they have different physical significances in each equation. 
When used to fit parts of the meson spectrum, the apparent success of the NRQM from this 
point of view is then due to its incorporation of variables numerically indistinguishable from 
their covariant versions together with a potential that fortuitously coincides (for a limited 
range of states) with a covariant one modified by an energy dependent constant term that 
varies from state to state. 

APPENDIX A: PAULI-FORM OF THE TWO-BODY DIRAC EQUATIONS 
FOR ct)+ = i)i+ AND THEIR RADIAL FORMS 

Reference [53] sets out Two-Body Dirac Equations containing general covariant interac- 
tions along with their accompanying Schrodinger-like forms. The general interactions consist 
of the eight Lorentz invariant forms corresponding to scalar, time and space- like vector stud- 
ied here along with five others: pseudoscalar, time- and space-like pseudovector, axial and 
polar tensor. When Eq.(2.26) is written in terms of the four four-component spinors ■01...4 
it decomposes into eight coupled equations. In [53] Long and Crater showed how these 
may be rearranged in Pauli-form or Schrodinger-like equations in terms of the combination 
0+ = V'l + i-^i ill the process providing a simpler coupling scheme than that used in [25] 
which involves coupled equations between ipi and ip^. Eq.(4.24) of reference [53] yields the 
following equation (simplified here for electromagnetic-like interactions {dJ = — —dG) 
and scalar interactions alone): 



+ 



1 



{M2Dt+ - MiD++) 



E1M2 + E2M1 



1 



[E^Dt^ + E^Dt^)]ct>. 



E1M2 + E2M1 



'+ 



{El - M2)0. 



'+ 



(A.l) 



in which the kinetic-recoil terms appear through the combinations: 




exp^ ai ■ p+ -a2 ■ d[L + g{l - ai ■ a2)] 
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expQ 



a2-p + -ai- d[L + g{l - ai ■ 0-2)] 



L>r+ = exp Q\ai-p + -a2-d [-L + g(l - ai ■ 02)] 



-Df" = exp^ CTi •p + -(72 ■ d[L -g{l + ai- (72)] 



(A.3) 



Manipulations using both sets of Pauli-matrices then lead to the form presented in the text 
in Eq.(2.35). 

We obtain the radial forms of Eq.(2.35) that wc use for our numerical solution for the 
general fermion-antifermion system by forming standard matrix elements of spin-dependent 
operators (see Appendix C of Ref.([25] )). We start from the general wave function of the 
form 

tpijm = ^ CiisRilsjyisjm] i = 1, 2, 3, 4 (A.4) 

l,s 

in which Rusj — is the associated radial wave function and yisjm is the total angular 
momentum eigenf unction. In terms of "D = E1M2 + E2M1 the corresponding radial forms 

then become 

s = 0, j^l 



{--^,+^^^^4^ + '2m^S + S' + 2e2A-A' 
dr r 

-{2g-iogiv) + g + Ly{-^-l) 

i V^(L + Ag) - 1(-L - 2^? + 2 log(P))'(-L - Agyjujoj 



, ^ . w(mi — mo) , „ 
+exp(-^? - L) ^ ^\ -g + L) 



,V j{j + 1) ,2/ N /A r\ 

— -^jij = b [wjujoj, (A.5) 



coupled to 



s = 1, j 



-(^-L-logp))';^-- 

dr r 



+\V^L +1(2 log(P) + (-L + 2g))'L']u^,^ 
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+exp(-^ - J) ^ {-y + L) — Ujoj = b [wjUjij, (A. 10) 

and s — — l + l 



2 



+ '^^j^^Y^ {2[g + ^]'|: + [(-^ - - 2j) + 3g]'l 

+V\L) - L'{log{V) - 2gy}u,+nj = b\w)u,_nj, (A.ll) 
coupled tos — — I — 1 

+[log(P)-2^-^(G + L)r^ 
[{3 + 1) log(P) - 2^((4i^ + 7j +A)g-g- L)]'\ 
+ + ^)'') - 27TT(^^^'^ + ^'(^^ - ^ - L - C)' + 1 log'(D)L')}«,+n, 

+ ^2i+Y^ ^^^^ + ^]'|: + [{-g - L){2j + 3) + 36;]'^ 

+2V'L + L'(log(P) - 2e?)'}ii,_n,- = b''{w)uj+nj- (A.17) 



APPENDIX B: NUMERICAL CONSTRUCTION OF MESON WAVE FUNC- 
TIONS 

We obtain from our computer program a numerical wave function -^(0;) normalized so 
that 

/+00 
u(xfdx = 1. (B.l) 
■00 
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The radial variable is related to a; by r = r^e^ and the radial wave function u{r) — 
■u(x)e~^/^/^/r^. Hence 



u{rYdr= / u{xYdx. (B.2) 

J-oo 

Now let Vn{r) be some radial basis functions that are orthonormalized so that 

r+oo 

/ Vn{r)Vn'{r)dr = 5nn'- (B.3) 

Thus 

oo 

u{r) = y^^UnVnjr) (B.4) 

n=0 

where 

/•+00 r+oo 

Un = Vn{r)u{r)dr = I Vn{x)u{x)dx. (B.5) 

Jq J-oo 

Note that Vn{x) = Vn{r)e^l'^ ^[tq so that we can compute the in a straightforward way. 
Thus we have as an approximation 

/' + OC 

''^('^) = "^nl^) / Vn{x)u{x)dx 

= ^ CnVn W = Wiv(r). (B.6) 

n=0 

Now we use a least squares fit to determine the c„ . In the limit of large N we have 
Cn — tin since we minimize the quantity 

/ + CX) 
\u{x)-WN{x)fdx (B.7) 
-oo 

For the VrSj^ we use harmonic oscillator (Laguerre) functions defined by 

vl{y)^c{n,k)e-y'l^y^Ll-'l\y') (B.8) 



in which c(n, k) = a{n+k-i/2)\ normalization constant and in terms oi z = y"^ 



„2„-fe+l/2 Jn 



So for example 
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Lr^\z) = 1 

Lf-'/2(^) = A; + l/2-z 

L^-'^'iz) = i[(5/2 + k - z)L\-''\z) - (1/2 + k)Ll-"\z) 
= [{k + 3/2)(A; + 1/2) - 2{k + 3/2)^ + ^^1/2 

Lt;f (^) = ^[(2n + 1/2 + A; - ^)L^V2(^) -{n + k- l/2)Ltf (^)] (B.IO) 
n + 1 

Thus letting y = r/a = ae^ we obtain 



vq{x) = c{0,k)a''exp{x{2k + l)/2) expi-a^e"^"" /2) 



'^'^''^ = ^lkTl/2'''^''^^^ + ~ "'"'"^^ 



""^^"^^ = + 1/2X^ + 3/2) ^°^"^^^^^ + ^/^^^^ + ^/^^ " + 3/2)a2e2- + aVl/2. 



^ \l {k + l/2)..{k + n-l/2) ""^^ ' {n -m)\{k- 1/2 + m)\m\^ ' 

(B.ll) 



APPENDIX C: TABLE VI - COMPARISON OF IMPORTANT FEATURES OF 
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APPROACHES TREATED IN THIS PAPER 



Wave Eqn 
Covariance 
Nonperturb. Tests 
# of Parametric fns 

Locality 



HC-PVA Wisconsin Iowa State Brayshaw Godfrey,Isgur 

Two-Body Dirac Reduced BSE Quasipotential Breit None 



Explicit Implicit 
Strng. ptnl -QED Wk ptnl. 
2 2 
101 5169 vs 73 

Local Non-local 



Running coupling cnst. Yes 



Yes 



Implicit Implicit Implicit 

Str. ptnl. Str. ptnl Str. ptnl. 

2 3 6 

RMS 50 vs 53 204 vs 111 85 vs 105 

Non-local Local Non-local 

Yes No Yes 
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TABLE I - 
NAMICS 

NAME 



MESON MASSES FROM COVARIANT CONSTRAINT DY- 

EXP. THEORY 



T : 


bb l^Si 


9.460( 


0.2) 


9.453( 


0.6 


T : 


bb I^Pq 


9.860( 


1.3) 


9.842( 


1.4 


T : 


bb l^Pi 


9.892( 


0.7) 


9.889( 


0.1 


T : 


bb l^Po 


9.913( 


0.6) 


9.921( 


0.5 


T : 




10.023( 


0.3) 


10.022( 


0.0 


T : 


bb 2^ Pq 


10.232( 


0.6) 


10.227( 


0.2 


T : 


bb 2^ Pi 


10.255( 


0.5) 


10.257( 


0.0 


T : 


bb 2^ P2 


10.269( 


0.4) 


10.277( 


0.8 


T : 


bb S^Si 


10.355( 


0.5) 


10.359( 


0.1 


T : 


bb A^Si 


10.580( 


3.5) 


10.614( 


0.9 


T : 


bb 5^Si 


10.865( 


8.0) 


10.826( 


0.2 


T : 


bb 6^Si 


11.019( 


8.0) 


11.013( 


0.0 


B : 


bu l^S'n 


5.279( 


1.8) 


5.273( 


0.1 


B : 


bd l^Sa 


5.279( 


1.8) 


5.274( 


0.1 


B* 


: bu l^Si 


5.325( 


1.8) 


5.321( 


0.1 


B. 


: bs l^So 


5.369( 


2.0) 


5.368( 


0.0 


B. 


: 6s l^S"! 


5.416( 


3.3) 


5.427( 


0.1 


'Ic 


cc l^So 


2.980( 


2.1) 


2.978( 


0.0 


lb ■ 


cc l^Si 


3.097( 


0.0) 


3.129( 


12.6 


Yn 


: cc l^Pi 


3.526( 


0.2) 


3.520( 


0.4 


Xo 


. 1 3 p 

. cc 1 


3.415( 


1.0) 




n A 


Xi 


: cc l^Pi 


3.510( 


0.1) 


3.507( 


0.2 


X2 


: cc 1^P2 


3.556( 


0.1) 


3.549( 


0.6 


Vc 


cc 2^So 


3.594( 


5.0) 


3.610( 


0.1 


i,: 


cc 2^ Si 


3.686( 


0.1) 


3.688( 


0.1 


^: 


cc l^Di 


3.770( 


2.5) 


3.808( 


2.0 


i^: 


cc 3^Si 


4.040( 10.0) 


4.081( 


0.2 




cc 2^Di 


4.159( 20.0) 


4.157( 


0.0 
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lb ■ 

Y ■ 


cc 3^ Di 


4.415( 6.0) 


4.454( 


0.4) 


D : 


CU l^Sn 


1 8651' 51 


1.866( 


0.0) 


D : 


Cd l^Sn 


1 869 f 5) 


1.873( 


0.1) 


D* 


: CU l^Si 


2 007f 5) 


2.000( 


0.4) 


D* 


: cd l^Si 


2 OlOf 5) 


2.005( 


0.3) 


D* 


: CU l^Pi 


2.422( 1.8) 


2.407( 


0.6) 


D* 


: cd 1^ Pi 


2.428( 1.8) 


2.411( 


0.5) 


D* 


: CU 1^P2 


2 4591' 2 0) 


2.382( 


11.3) 


D* 


: cd 1^ P2 


2 459f 4 0) 


2.386( 


3.5) 


-Do 


: cs l^So 


1 968f 6) 


1.976( 


0.5) 


D* 

^ s 


: cs l^Si 


2.112( 0.7) 


2.123( 


0.9) 


D* 


: cs l^Pi 


2 535f 3) 


2.511( 


6.2) 


D* 


■ cs 1^ Po 


2.574( 1.7) 


2.514( 


9.6) 


K : 


SU l^Sn 


4941' 01 


0.492( 


0.0) 


K : 


Sd l^Sn 


498f 01 


0.492( 


0.4) 


K* 


: SU l^Si 


892f 21 


0.910( 


0.6) 


K* 


■ sd l^S, 


896f 31 

vy • vy »_y vy \ vy • vy / 


0.910( 


0.3) 


-•^ 1 


: SU l^Pi 


1.273( 7.0) 


1.408( 


3.2) 




: SU 1^ Pq 


1.429( 4.0) 


1.314( 


0.7) 


K^ 


: SU l^Pi 


1.402( 7.0) 


1.506( 


1.0) 




: SU 1^ P2 


1.425( 1.3) 


1.394( 


0.5) 




: sd 1^ P2 


1.432( 1.3) 


1.394( 


0.6) 


K* 


■ SU 2^ Sn 

• ^ K-^ I J 


1 460f 30 01 

_i_ • J- vy vy \ vy vy • vy / 


1.591( 


0.2) 


K* 


: SU 2^ Si 


1.412( 12.0) 


1.800( 


6.7) 




: sir 1^D2 


1 773( 8 01 

-1- ■ 1 1 vy I vy ■ vy j 


1.877( 


0.8) 


K* 


: SU 1^ 


1 714(' 20 01 


1.985( 


1.4) 




: SU 1^D2 


1 816f 10 01 

J- • vy J- vy I J. vy • vy 1 


1.945( 


1.3) 




: SU 1^ Ds 


1 770f 10 01 

ivy\ xvy«vy/ 


1.768( 


0.0) 




. SU Oq 


i.OOUl^ OU.U j 


z.ioo 


1 A\ 
i 1 


k; 


: sU 23P2 


1.975( 22.0) 


2.098( 


0.2) 


Kl 


: SU 1^F4 


2.045( 9.0) 


2.078( 


0.1) 


K2 


: SU 2^D2 


2.247( 17.0) 


2.373( 


0.5) 
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■ SU l^Gr. 


2 382( 33 01 


2.344( 


0.0) 


3 


: sTI 2^ F3 


2.324( 24.0) 


2.636( 


1.9) 




: 2^ F4 


2.490( 20.0) 


2.757( 


1.6) 


6 • 


ss l^Si 


1.019( 0.0) 


1.033( 


2.2) 


fn 


ss 1^ Pq 


1.370( 40.0) 


1.319( 


0.0) 




: ss 1^ Pi 


1.512( 4.0) 


1.533( 


0.3) 




: ss 1^P2 


1.525( 5.0) 


1.493( 


0.3) 




ss 2^Si 


1.680( 20.0) 


1.850( 


0.8) 


6 ■ 


ss l^Ds 


1 854f 7 0) 


1.848( 


0.0) 


J ^ 


ss 2^ P2 


2.011( 69.0) 


2.160( 


0.1) 


J -i 


ss 3^ P2 


2 297f 28 0) 


2.629( 


1.6) 


TT : 


ud I^jSq 


0.140( 0.0) 


0.144( 


0.2) 


o : 


ud l^Si 


0.767( 1.2) 


0.792( 


0.1) 


^1 


lid 1^ Pi 


1 2311' 10 01 


1.392( 


2.1) 


(In 


• ud 1^ Pn 


1 450f 40 Ol 


1.491( 


0.0) 


Oi 


: lid 1^ Pi 


1.230( 40.0) 


1.568( 


0.7) 




: lid 1^ P2 


1.318( 0.7) 


1.310( 


0.0) 


TT : 


ud 2^Sq 


1.300( 100.0) 


1.536( 


0.1) 


iO : 


ud 2^ Si 


1.465( 25.0) 


1.775( 


1.4) 


TTn 
2 


• ud l^Do 


1 670f 20 0) 


1.870( 


0.9) 


H ■ 


\JJ\AJ J- -I—-' 1 


1 700( 20 01 


1.986( 


1.9) 


r J 


: ud V'D-i 


1.691( 5.0) 


1.710( 


0.0) 


TT : 


ud 3^5*0 


1.795( 10.0) 


2.166( 


7.9) 




ud 3^5*1 


2.149( 17.0) 


2.333( 


0.7) 


Oa 


: ud 1^P4 


2.037( 26.0) 


2.033( 


0.0) 


'' z 


: ud 2^D2 


2.090( 29.0) 


2.367( 


0.5) 


P3 




Z.ZO\j[ 40. u ) 


2.305( 


0.0) 


P5 


: ud l^Gs 


2.330( 35.0) 


2.307( 


0.0) 


P6 


: l^He 


2.450( 130.0) 


2.547( 


0.0) 


x' 




0.0 


101.0 



TABLE II COMPARISON OF MESON MASSES FROM 

WISCONSIN MODEL II and COVARIANT CONSTRAINT DYNAMICS 

NAME EXP. WISC2 CCD 



T : bb l^Si 
T : bb l^Po 
T : bb l^Pi 
T : bb 1^P2 
T : bb 2^Si 
T : bb 2^Po 
T : bb 2^Pi 
T : bb 2^P2 
T : bb 3^Si 
B : bd I^Sq 
r]c : cc I'^So 
: cc l^^i 
Xo : cc l^Po 
Xi : cc l^Pi 
X2 : cc 1^P2 
?7c : cc 2^So 
V' : cc 2^,51 
D : cd l^,So 
D* : cd l^^i 
P's : cs I^Sq 
D* : cs l^Si 
K : sd I^Sq 
K* : sd l^Si 
(j) : ss l^Sx 
: ss 23^1 



Q 4601' 

U .t:\JW 1 


2\ 


.OWL/ 1 


1 Si 


Q 8Q2(' 


n 71 

W. f j 




n 61 


1 n n2S(' 


Si 


1 D 9'^9( 


1 11 


1 n 9^^( 


fil 

W . W 1 


1 n 2fi8f 


n 61 

W.W 1 


in 3551' 


n 51 


5 27Q( 


2 11 


2 Q7Qf 


1 91 




11 


s 41 


1 01 

X . W 1 




11 

W. X ( 


3 556 f 

kJ •\J\J\J I 


11 

W. X 1 


5Q4(' 


5 nl 

W.W j 


686/' 

(J .WOW I 


n 11 

W. X 1 


1.869( 


0.5) 


2.010( 


0.6) 


1.969( 


0.7) 


2.110( 


2.0) 


0.498( 


0.0) 


0.896( 


0.3) 


1.019( 


0.0) 


1.680( 50.0) 



0.0 



9.426( 62.6) 
9.862( 0.1) 
9.892( 0.0) 
9.917( 0.7) 
10.028( 1.3) 
10.238( 1.5) 
10.256( 0.0) 
10.270( 0.2) 
10.359( 0.7) 
5.381( 137.2) 
2.967( 1.4) 
3.167( 272.4) 
3.402( 5.1) 
3.493( 17.5) 
3.548( 4.0) 
3.621( 1.5) 
3.668( 17.9) 
1.983( 574.6) 
2.010( 0.0) 
2.097( 671.1) 
2.148( 52.7) 
0.743(3340.4) 
0.870( 5.1) 
1.019( 0.0) 
1.510( 0.9) 
5168.9 



9 454 


' 2 01 


9 845 


' 4 51 


9 890 


' ll 

W . X 1 


9 91 9( 


' 1 61 

X . W J 


10.024 


' 0.1) 


1 229 


' 41 


1 257 


' 91 


10.276 


[ 3.1) 


1 359 

X W .(JWC 


' 71 

W.I \ 


5 274( 


' 31 

W.W \ 


2 975 


'oil 

W . X 1 


3 1201' 

W . XZj W 1 


28 81 


3 41 2 

W . ^ X ^ 


' 21 

W . ^ 1 


3 505 

(J. WW w 


' 1 81 

X . O 1 


3 538 f 

W. WWO I 


18 11 

XO . X 1 


3 61 1 

W. WX X 


' 61 

W.W \ 


3 6881 

W. WOO 


' 31 

W.W I 


1 875 

X . O 1 w 


' 1 51 

X . W ( 


2.003 


' 1-9) 


1.968 


; 0.1) 


2.106 


; 0.6) 


0.498( 


; 0.0) 


0.918 


; 3.5) 


1.020 


; 0.0) 


1.424 


; 2.1) 




72.8 
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TABLE III -COMPARISON OF MESON MASSES FROM 
SPENCE-VARY MODEL and COVARIANT CONSTRAINT DYNAMICS 

NAME EXP. SPENCE & VARY CCD 



T : 


bb l^Si 


9.460(9.460) 


9.452(-8) 


9.444(-16 


T : 


bb l^Po 


9.860(9.860) 


9.843(-17) 


9.836(-24 


T : 


bb l^Pi 


9.892(9.893) 


9.863(-29) 


9.886(-7 


T : 


bb 1^P2 


9.913(9.913) 


9.928(+15) 


9.921(+8 


T : 


bb 2^Si 


10.023(10.023) 


9.996(-27) 


10.022(+1 


T : 


bb 2^Po 


10.232(10.232) 


10.198(-34) 


10.230(+2 


T : 


bb 2^ Pi 


10.255(10.255) 


10.214(-41) 


10.261(+6 


T : 


bb 2^P2 


10.268(19.269) 


10.270(+2) 


10.284(+17 


T : 


bb ?,^S-i_ 


10.355(10.355) 


10.331(-24) 


10.367(+12 


T : 


bb A^Si 


10.580(10.580) 


10.611(+31) 


10.627(+47 


T : 


bb b'^Si 


10.865(10.865) 


10.860(-5) 


10.645(-20 


T : 


bb G^Si 


11.019(11.019) 


11.086(+67) 


11.036(17 


B : 


bu I^Sq 


5.271(5.279) 


5.342(+63) 


5.267(-12 


B* 


: bu l^Si 


5.352(5.325) 


5.347(-5) 


5.317(-8 


IC 


cc V-Sq 


2.979(2.980) 


2.993(+14) 


2.969(-ll 


lb : 

T • 


cc l^Si 


3.097(3.097) 


3.091(-6) 


3.128(+31 


Yn 


: cc l^Pi 


3.526(3.526) 


3.471(-55) 


3.520(-6 


Yn 


: cc l^Po 


3.415(3.415) 


3.383(-32) 


3.396(-19 


Yi 
Ai 


: cc l^Pi 


3.511(3.511) 


3.461(-50) 


3.504(-7 


Y9 


: cc 1^P2 


3.556(3.556) 


3.556(0) 


3.555(-l 




CC Z Oq 






O.DUD^^-^iZ 


i;: 


CC 2^Si 


3.686(3.686) 


3.688(+2) 


3.688(+2 


i;: 


cc l^Di 


3.770(3.770) 


3.741(-29) 


3.806(+36 


V^: 


cc 3^Si 


4.040(4.040) 


4.104(+64) 


4.083(+43 


^: 


cc 2^Di 


4.159(4.159) 


4.136(-23) 


4.161(+2 


^: 


cc 3^Di 


4.415(4.415) 


4.456(+41) 


4.462(+47 


D 


cu I^Sq 


1.865(1.8645) 


1.897(+32) 


1.854(-10 


D* 


: cu l^Si 


2.007(2.007) 


2.004(-3) 


1.991(-16 



69 



D* : cu l^Pi 
Ds : cs I^Sq 
Dl : cs l^Si 
K : su I^Sq 
K* : su l^Si 
Ki : su l^Pi 
: su l^Pi 
^2* : su 1^P2 
K2 : 1^D2 
(p : l^^i 

/2 : 13P2 

: ss 2^Si 
TT : ud l^So 
p : ud l^Si 
bi : ud l^Pi 

01 : ud l^Pi 

02 : ud 1^P2 
TT : ud 2^So 

772 : liC^ 1^-02 

RMS 



2.420(2.422) 
1.971(1.969) 
2.110(2.112) 
0.494(0.494) 
0.892(0.892) 
1.270(1.273) 
1.406(1.402) 
1.426(1.426) 
1.770(1.776) 
1.019(1.019) 
1.525(1.525) 
1.680(1.680) 
0.140(0.140) 
0.768(0.769) 
1.232(1.230) 
1.260(1.230) 
1.318(1.318) 
1.300(1.300) 
1.670(1.670) 
0.0 



2.358(-72) 
1.968(-3) 
2.076(-34) 
0.495(+l) 

0.916(+24) 

1.287(+17) 
1.330(-76) 
1.330(-96) 

1.633(-137) 
1.020(+1) 
1.526(+1) 
1.645(-35) 
0.135(-5) 

0.812(+44) 
1.219(-13) 
1.223(-37) 

1.367(+49) 
1.439(+139) 

1.515(-155) 
50 



2.373(-47) 
1.981(+12) 
2.137(+25) 
0.511(+17) 
0.887(-5) 
1.327(+57) 
1.405 (+3) 
1.348(-78) 
1.709(-85) 
1.048(+29) 
1.488(-37) 
1.803(+123) 
0.143(+3) 
0.736(-33) 
1.255(+25) 
1.534(+185) 
1.223(-95) 
1.474(174) 
1.780(+110) 
53 
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TABLE IV - COMPARISON OF MESON MASSES FROM 
BRAYSHAW MODEL and COVARIANT CONSTRAINT DYNAMICS 

NAME EXP. BRAYSHAW CCD 



T : 


bb l^Si 


9.460( 


0.2) 


9 452f 1 3) 


9.451( 


1.7 


T : 


bb I^Pq 


9.860( 


1.3) 


9 866f 3) 


9.842( 


2.5 


T : 


bb l^Pi 


9.892( 


0.7) 


9 910( 4 51 


9.889( 


0.1 


T : 


bb l^Po 


9.913( 


0.6) 


9 926( 2 5) 


9.920( 


0.7 


T : 


bb 2^ Si 


10.023( 


0.3) 


10 007f 4 8) 


10.023( 


0.0 


T : 


bb 2^Pq 


10.232( 


0.6) 


10.214( 4.9) 


10.229( 


0.1 


T : 


bb 2^ Pi 


10.255( 


0.5) 


10.252( 0.1) 


10.258( 


0.1 


T : 


bb 2^P2 


10.268( 


0.4) 


10.265( 0.2) 


10.278( 


1.8 


T : 


bb S^Si 


10.355( 


0.5) 


10.342( 2.8) 


10.360( 


0.4 


T : 


bb A^Si 


10.580( 


3.5) 


10 6621' 9 41 


10.617( 


1.9 


B : 


bu l^Sn 


5.279( 


1.8) 


5 332 f 13 7) 


5.270( 


0.3 


B* 


: bu \^Si 


5.325( 


1.8) 


5 377f 13 2) 


5.317( 


0.3 


'Ic ■ 


cc V-Sq 


2.980( 


2.1) 


3 Ollf 3 51 


2.976( 


0.0 


lb ■ 


cc l^Si 


3.097( 


0.1) 


3 129f 21 01 


3.127( 


17.8 


Yn 


: cc l^Pi 


3.524( 


0.2) 


3.498( 13.0) 


3.520( 


0.3 


Yn 


: cc l^Po 


3.415( 


1.0) 


3.410( 0.3) 


3.409( 


0.4 


Yi 
Al 


: cc 1^ Pi 


3.510( 


' 0.1) 


3 514(' 21 


3.508( 


0.2 


Yo 


■ cc 1^ P9 


3.556( 


0.1) 


3 540f 5 21 


3.547( 


1.5 


'Ic ■ 


cc 2^Sr, 

L- ^ LJ (J 


3.594( 


5.0) 


3 580f 21 


3.612( 


0.3 


lb : 


CC 2^ Si 


3.686( 


0.1) 


3.680( 0.7) 


3.691( 


0.4 


ih ■ 

w ■ 


CC i Ui 


770/ 


Z.O j 


0.1 i oy u.u J 


3.811( 


4.0 




CC 2>^Si 


4.040( 


10.0) 


4.246( 8.0) 


4.086( 


0.4 




cc 2^Di 


4.159( 


20.0) 


4.288( 0.8) 


4.163( 


0.0 


D : 


cu I^Sq 


1.865( 


0.5) 


1.903( 24.2) 


1.864( 


0.0 


D* 


: cu l^Si 


2.007( 


1.4) 


2.046( 24.5) 


1.997( 


1.7 


D* 


: cu l^Pi 


2.422( 


1.8) 


2.428( 0.1) 


2.413( 


0.3 


D* 


: cu l^Pa 


2.459( 


2.0) 


2.458( 0.0) 


2.383( 18.8 


Ds 


: cs I^Sq 


1.969( 


0.6) 


1.976( 0.8) 


1.974( 0.4 



71 



D* 

^ s 


: cs l^Si 


2.112( 2.0) 


D* 

s 


: cs l^Pi 


2.535( 0.3) 


D* 

^ s 


: cs 1^ P2 


2.574( 1.7) 


K 


su ^^Sq 


0.494( 0.0) 


K* 


: su l^Si 


0.892( 0.2) 


i^i 

^ 1 


: su l^Pi 


1.273( 7.0) 




: su 1^ Po 


1.430( 4.0) 


i 


: siT 1^ Pi 


1.402( 7.0) 




: su 1^P2 


1.425( 1.3) 


K* 


: su 1^ Dl 


1.714( 20.0) 


^ z 


: su 1^ 


1.816( 10.0) 


^ o 


: l^Ds 


1.770( 10.0) 


6 : 


ss l^Si 


1.019( 0.0) 


JO ■ 


ss I^Pq 


1 370C 40 Q) 


fi : 


ss l^Pi 


1.512( 4.0) 




ss 1^P2 


1.525( 5.0) 


6 : 


SS 2^ Si 


1.680( 20.0) 




ss l^D:i 


1.854( 7.0) 


TT : 


ud V-Sq 


0.140( 0.0) 


O ' 
r • 


ud l^Si 


0.767( 1.2) 


h^ ■ 
w\ . 


ud 1^ Pi 


1 2311' 10 0) 




ud l^Po 


1.450( 40.0) 




ud l^Pi 


1.230( 40.0) 


do 


ud 1^ P2 


1.318( 7.0) 


TT : 


ud 2^So 


1.300( 100.0) 




ua i U2 


1 dTCif on n^ 

i.D / U\^ ZV.y) ) 


p ■ 


ud l^Di 


1.700( 20.0) 


P3 


ud l^Da 


1.691( 5.0) 


x' 




0.0 



2.134( 


6.6) 


2.119 


[ 0.7) 


2.515( 


7.2) 


2.515 


' 7-0) 




fi^ 


2 51 3(' 


17 01 


0.495( 


0.0) 


0.492 


' 0.1) 


905 ( 


51 


908 


' 71 


1.355( 


1.1) 


1.421 


' 3.6) 


1.086( 


10.8) 


1.349 


' 0.6) 


1.294i 


3.4) 


1.524 


' 4-3) 


1 4091 


21 


1 399 


' 51 


1 fi90( 


01 


2 004 


' 2 61 


1.764( 


0.4) 


1.892 


' 0.8) 


1.770( 


0.0) 


1.780 


' 0.0) 


1.022( 


0.1) 


1.030 


\ 2.1) 


1 185( 


41 


1 345 


' 01 


1.446i 


4 51 


1 546 


' 1.2) 


1.511 


1 1 


1 496 


' 41 


1.778( 


0.4) 


1.860 


[ 1.4) 


1.922( 


1.4) 


1.856 


' 0.0) 


0.140( 


0.0) 


0.143 


' 0.2) 


776 ( 


01 


790 


' 21 


1.202i 


11 


1.411 


' 4.4) 


990( 


2.4) 


1.542 


' 11 


1.253( 


0.0) 


1.590 


[ 1.3) 


1.302( 


0.2) 


1.318 


' 0-0) 


1.028( 


0.1) 


1.543 


' 0.1) 


1.593( 


0.2) 


1.883 


[ 1.6) 


1.741( 


0.1) 


1.998 


[ 3.4) 


1.680( 


0.0) 


1.722 


; 0.2) 


204.2 




111.0 
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TABLE V - COMPARISON OF MESON MASSES FROM 
ISGUR-WISE MODEL and COVARIANT CONSTRAINT DYNAMICS 

NAME EXP. ISGUR&WISE CCD 



T : 


bb l^Si 


9.460( 


0.2) 


9.460( 


0.0) 


9.453( 


0.8 


T : 


bb I^Pq 


9.860( 


1.3) 


9.850( 


0.5) 


9.842( 


1.6 


T : 


bb l^Pi 


9.892( 


0.7) 


9.880( 


1.4) 


9.889( 


0.1 


T : 


bb l^Po 


9.913( 


0.6) 


9.900( 


1.8) 


9.921( 


0.6 


T : 


bb 2^ Si 


10.023( 


0.3) 


10.000( 


6.9) 


10.023( 


0.0 


T : 


bb 2^ Po 


10.232( 


0.6) 


10.230( 


0.0) 


10.228( 


0.2 


T : 


bb 2^ Pi 


10.255( 


0.5) 


10.250( 


0.3) 


10.257( 


0.0 


T : 


bb 2^ P2 


10.269( 


0.4) 


10.260( 


1.0) 


10.277( 


0.8 


T : 


bb ?>^Si 


10.355( 


0.5) 


10.350( 


0.3) 


10.359( 


0.2 


T : 


bb A^Si 


10.580( 


3.5) 


10.630( 


2.4) 


10.615( 


1.2 


T : 


bb 5^Si 


10.865( 


8.0) 


10.880( 


0.0) 


10.828( 


0.2 


T : 


bb 6^Si 


11.019( 


8.0) 


ll.lOOf 


1.2) 


11.014( 


0.0 


B : 


bv l^S'n 


5.279( 


1.8) 


5.310( 


3.3) 


5.272( 


0.2 


B* 


: bu l^Si 


5.325( 


1.8) 


5.370( 


6.9) 


5.319( 


0.1 


B. 


: bs l^So 


5.369( 


2.0) 


5.390( 


1.2) 


5.368( 


0.0 


B. 


: bs l^Si 


5.416( 


3.3) 


5.450( 


1.4) 


5.426( 


0.1 


'Ic 


cc l^So 


2.980( 


2.1) 


2.970( 


0.2) 


2.978( 


0.0 


lb ■ 


cc l^Si 


3.097( 


0.0) 


3.100( 


0.1) 


3.128( 


14.1 


Vn 


: cc l^Pi 


3.526( 


0.2) 


3.520( 


0.5) 


3.520( 


0.5 


Vn 


: cc 1^ Po 


3.415( 


1.0) 


3.440( 


4.4) 


3.408( 


0.4 


Xi 


. ^7; 1 3 p 
. CC 1 ri 


3.510( 


0.1) 


O.OIUI 


U.U ) 


3.507( 


0.2 


X2 


: cc 1^P2 


3.556( 


0.1) 


3.550( 


0.5) 


3.548( 


0.9 


Vc 


cc 2^So 


3.594( 


5.0) 


3.620( 


0.4) 


3.611( 


0.2 


V^: 


cc 2^Si 


3.686( 


0.1) 


3.680( 


0.5) 


3.689( 


0.1 


^: 


cc l^Di 


3.770( 


2.5) 


3.820( 


4.2) 


3.809( 


2.5 


^: 


cc 3^ Si 


4.040( 10.0) 


4.100( 


0.5) 


4.082( 


0.2 


i^: 


cc 2^Di 


4.159( 20.0) 


4.190( 


0.0) 


4.159( 


0.0 




cc 3^Di 


4.415(6.0) 


4.450( 


0.4) 


4.456(0.6 
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D : 


CU l^Sn 

-M- Aw/ II 


1 865( 5) 


1.880( 2.7) 


1 865( 0) 


D* 


: CU l^Si 


2 007( 5) 


2.040( 12.6) 


1 998( 8) 

_i_ ■ t_/ »_/ yy \ vy ■ yy 1 


D* 


: CU l^Pi 


2 422f 1 8) 


2.440( 


0.9) 


2 408f 6) 


D* 


: CU l^P2 


2.459( 2.0) 


2.500( 


3.8) 


2 381( 13 6) 

• ■ y yy \ ■ y • yj 1 


S 


: cs l^So 


1 968f 6) 


1.980( 


1.4) 


1 976f 6) 

• t_y 1 yj \ yj * yj 1 


D* 


: cs l^Si 


2.112( 0.7) 


2.130( 


3.0) 


2.121( 0.8) 


D* 


: cs 1^ Pi 


2.535( 0.3) 


2.530( 


0.4) 


2.512( 6.7) 


D* 


■ cs l^Po 


2.574( 1.7) 


2.590( 


0.9) 


2 513(11 6) 


K : 


SU l^Sn 


494f 0) 


0.470( 


8.0) 


494f 0) 

Vy • ^ JI \ Vy • Vy / 


K* 


: SU l^Si 


892f 2) 


0.900( 


0.1) 


907f 5) 

vy • yj 1 \ Vy • Vy / 


*■ i 


: SU l^Pi 


1.273( 7.0) 


1.340( 


0.5) 


1.411( 2.2) 





: SU I^Pq 


1.429( 4.0) 


1.240( 


2.3) 


1.323( 0.7) 
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